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Abstract 

Future detection of Extensive Air Showers (EAS) produced by Ultra High Energy Cosmic Particles by means of space 
based fluorescence telescopes will open a new window on the universe and allow cosmic ray and neutrino astronomy 
at a level that is virtually impossible for ground based detectors. In the context of the Extreme Universe Space 
Observatory (EUSO) project, an end-to-end simulation of EAS observation with a spatial detector has been designed 
(EUSO Simulation and Analysis Framework, ESAF). This paper describes the detailed Monte-Carlo developed to 
simulate all the physical processes involved in the fluorescence detection technique, from the EAS development to 
the instrument response. Particular emphasis is given to modeling the light propagation in the atmosphere and the 
effect of clouds. The simulation is used to assess the performances of EAS spatial detection. Main results on energy 
threshold and resolution, direction resolution and X max determination are reported. Results are based on EUSO 
telescope design, but are also extended to larger and more sensitive detectors. 

Key words: cosmic rays; space detectors; extensive air showers; atmosphere models; 
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1. Introduction 

Ultra-High Energy Cosmic Particles (UHECP) 
with energies in excess of ~ 10 19 eV hit the Earth 
with a very low flux which falls down to less of 
one particle km -2 sr _1 century -1 for particles with 
energy E > 10 20 eV [1]. 
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The observation of UHECP and the interpreta- 
tion of the related phenomenology is one of the most 
challenging topics of modern High-Energy Astro- 
Particle Physics. Direct detection is impossible at 
these energies, due to the exceedingly low flux, but 
UHECP can be detected by observing the Exten- 
sive Air Showers (EAS) produced by the interaction 
of the primary particle with the Earth atmosphere 
(see [2,3,4], and references therein, for recent reviews 
on these topics). 

The ground-based Pierre Auger Observatory 
(PAO) [5,6,7,8] is currently taking data: its south 
site (surface of 3000 km 2 ) in Argentina and its 
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forthcoming north site in the US (20000 km 2 ) will 
provide, in the next few years, a clear understanding 
of many important topics [1,9,10,11,12]. However, a 
next generation of experiments with sensitivity to 
UHECP flux well above the PAO one at the highest 
energies is essential but it seems to be difficult to 
realize by extending the ground-based technique at 
ever larger arrays. 

Ground-based experiments detect EAS either by 
sampling charged particles at ground by means of 
suitable counters and/or by detecting the fluores- 
cence light emitted by the atmospheric nitrogen 
when excited by the charged particles in the shower. 
While the former technique is intrinsically ground 
based, the latter is not, being the fluorescence light 
emission isotropic. However, fluorescence light de- 
tection is limited to moonless nights and is affected 
by weather conditions. 

In 1979 J. Linsley [13] proposed for the first time 
the idea of detecting EAS fluorescence light with 
a telescope in orbit around the Earth. While this 
approach is technically, programmatically and eco- 
nomically challenging, it is one of the most promis- 
ing that might allow a big step forward in this field, 
opening a new window on the universe. A space 
based experiment might observe EAS on an effective 
surface of a few times 10 5 km 2 , increasing the event 
statistics by more than an order of magnitude com- 
pared to existing and running ground based experi- 
ments, even when taking into account that the duty 
cycle for observing EAS with fluorescence is limited 
to clear and moonless nights. Moreover, the huge 
target mass of the observed air volume (> 10 12 kg), 
three orders of magnitude more than the current 
or foreseen projects in the South Pole ice [14] or 
mediterranean sea [15] ) may open the era of high 
energy neutrino astronomy, at least above energies 
of the order of 10 19 eV or so. However, we do not 
discuss high energy neutrino detection in this paper. 

The relevance of this innovative approach is 
widely recognized by the scientific community. 
The Extreme Universe Space Observatory (EUSO) 
project [16], originally proposed in 2000 as response 
to F2/F3 European Space Agency (ESA) call for 
mission successfully completed the phase A concep- 
tual study. Nevertheless it did not continue only 
for programmatic reasons related to International 
Space Station and Space Shuttle availability. The 
inclusion of UHECP physics in the ESA Cosmic 
Vision 2015-2025 program [17] provides a suitable 
framework for the study of future space missions. 
An introduction to the basic issues concerning the 



detection of UHECP with a space detector can be 
found in [18]. 

In this paper we therefore report the results ob- 
tained in the context of the EUSO collaboration by 
means of a full simulation and reconstruction soft- 
ware code named ESAF, standing for EUSO Simu- 
lation and Analysis Framework. 

Particularly, we report about results concerning 
energy threshold as function of UHECP direction 
and detector size and geometry, direction and en- 
ergy resolution, the dependence on atmospheric con- 
ditions and clouds. 

This paper is structured as follows: section 2 
presents the ESAF design; section 3 describes the 
simulation part, giving details on each modules, 
from the EAS development to the detector re- 
sponse; section 4 contains a description of the re- 
construction tasks; section 5 presents several results 
obtained with the ESAF code; the most relevant 
features are summarized in the conclusions. 

2. The ESAF design 

ESAF handles the whole simulation and recon- 
struction chain from primary particle interaction in 
the Earth atmosphere until the final reconstruction 
of the event. The code includes: 

- extensive air shower simulation both by means 
of internally developed algorithms (e.g. Greizen- 
Ilina-Linsley (GIL) parametrization, [19]) and 
with interfaces to existing widely used codes (e.g. 
CORSIKA, [20]); 

- a complete description of the atmosphere, includ- 
ing aerosols 3 and clouds 

- fluorescence and Cherenkov light production 

- a complete simulation of photon propagation, 
from production point up to the telescope, in- 
cluding detailed interactions with ground and at- 
mosphere, and a Monte-Carlo code dealing with 
multiple scattering 

- telescope optics simulation 

- telescope geometry 

- telescope photodetector simulation 

- electronics and trigger simulation 

- background simulation 

- pattern recognition and shower signal identifica- 
tion above background 

- reconstruction of direction, energy, and the 
so called slant depth of the shower maximum 

3 Aerosols are airborne liquid droplets or solid particles, 
which may exist at low altitude. 
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(J max ), i.e. the distance (expressed in g/cm 2 of 
traversed material) of the shower maximum from 
the top of the atmosphere. 

Although the simulation was done having in mind 
the specific EUSO design, the flexibility of the ESAF 
code allowed us to obtain results that are quite gen- 
eral and might be useful for any future project. 




Figure 1. The upper level structure of ESAF. The simulation 
and reconstruction modules are detailed in Fig. 2 and Fig. 3 
respectively. 




EAS development in atmosphere Signal generation in the detector 




Figure 2. ESAF simulation structure. 



2.1. General structure of the ESAF code 

The ESAF code is written mostly in C++ with 
a few external Fortran libraries. A complete object- 
oriented design is used to allow flexibility, modular- 
ity and maintenability of the code. The code is di- 
vided in two independent sections, simulation and 
reconstruction, that share a common infrastructure 
and a set of common libraries. The code is not public, 
but is available upon request to interested users 4 . 

4 http:/ /www. ge.infn.it/~pesce/esaf.html 




Figure 3. ESAF reconstruction structure. 

The main structure of the program and the inter- 
actions between the two main sections are shown in 
Fig. 1. 

The simulation is structured in 6 main sections 
(see Fig. 2): shower simulation starting from the 
expected UHECP flux, fluorescence and Cherenkov 
photons production in atmosphere, light propaga- 
tion from the production point up to the orbiting 
telescope with simulation of occurring photon inter- 
actions during the propagation, simulation of tele- 
scope optics (lenses and focal surfaces), simulation 
of front end electronics and trigger. In this paper 
we do not cover implementation details that can be 
found in [21,22,23]. 

The output of the simulation module is a ROOT 
file [24] with a structure that contains the same in- 
formation that is expected from real data plus a set 
of "Monte Carlo truth" data. The level of detail is 
user- configurable, to optimize the output size to the 
goal of the simulation (i.e. physics, optics, trigger, 
etc.). To provide the user with an easy access to the 
data, a set of detector event viewers is provided with 
ESAF. They allow a coherent display of the run in- 
formations, histograms and several 2D and 3D views 
of the event development in the detector and on the 
focal surface (some example plots are in section 3.4). 

The reconstruction module reads back the ROOT 
file representing the fully simulated shower with 
background and reconstructs the events as if they 
were real data events from the telescope. Several 
modules allow pattern recognition and identifica- 
tion of EAS signal above incoherent star light and 
moonlight background, 3D direction reconstruction, 
energy and X max reconstruction. The reconstruc- 
tion structure is shown in Fig. 3. 

The reconstruction framework is the main struc- 
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ture that acquires the data (real or simulated) from 
an input module and builds the chain of modules 
necessary to reconstruct the events. The event re- 
construction is divided into different tasks (e.g. 
background subtraction, direction reconstruction, 
. . .) and for each task different modules may be 
available. The use of such a structure allows to eas- 
ily replace or exclude a given module. In this way 
different algorithms or algorithm combinations can 
be tested and compared on the same events. As 
in the simulation part, the reconstruction modules 
are configurable by the user. The event container 
stores all the relevant informations about the recon- 
structed event. Every module can access the event 
container, reads the existing data stored in it and 
writes its own results. 

2.2. Atmosphere description 

In the simulation process, a description of the at- 
mosphere is required for proper simulation of shower 
development, light production and photons transfer 
toward the telescope entrance pupil (see Fig. 2). The 
atmosphere parameters are also required during the 
reconstruction process. 

In ESAF, the atmosphere and Earth surface are 
spherically shaped, and several models have been in- 
troduced, to form database of atmospheric profiles: 
the US-Standard model and its complements [25] 
and the NRLMSIS-00 (MSISE) empirical model [26]. 
The formers provide temperature, pressure, refrac- 
tive index and density profiles, up to 100 km, as well 
as ozone and water vapour profiles. Available data 
are mean values over time of atmospheric data at 
different latitudes. With MSISE model, atmospheric 
profiles are provided for given latitude, longitude, 
day of the year and hour. 

Other crucial data concerning atmosphere de- 
scription have been implemented in ESAF. The 
presence of clouds and aerosols has an impact on 
the UHECP detection from space. The cloud and 
aerosol effects have to be evaluated with the simu- 
lation. The cloud mean fraction in the field of view 
of an orbiting telescope is around 70%. To describe 
cloud main characteristics, informations collected 
with the satellite TOVS (TIROS Operational Ver- 
tical Sonder) Path-B [27] are used. To describe 
aerosols ESAF is also supplied with several atmo- 
spheric models and parameters which exist in the 
L0WTRAN7 5 software [28]. 

5 L0WTRAN7 is a standard program used in the atmosphere 



A set of dedicated configuration files allow to 
choose one particular atmosphere model and sta- 
tus among all the possibilities available in ESAF, 
before runs. The atmospheric parameters are then 
used at each step where they are required during 
the simulation or the reconstruction process. 

3. The simulation code 

3.1. Shower simulation 

An accurate and physically reliable simulation of 
EAS is fundamental for the correct understanding of 
the UHECP phenomenology and the possible exper- 
imental outcomes. During the last 20 years several 
authors have dealt with the very demanding task 
of preparing detailed UHECP simulation programs. 
Particularly, CORSIKA [20] and AIRES [29] , are 
at present widely used and represent the 'state of 
the art' in atmospheric shower simulation. Another 
shower simulator, UNISIM [30], was developed by 
some ESAF developers, and has the unique feature 
of simulating neutrino interactions in atmosphere. 

The basic ESAF approach, therefore, does not 
uniquely consist in the development of brand new 
shower simulators, but also in having the capabil- 
ity to handle external shower simulation programs. 
This attitude is realized by implementing in ESAF a 
generic interface by means of a 'Shower Track' object 
which contains the structure for an extremely de- 
tailed description of a single shower event, a descrip- 
tion that is even more accurate than any of the out- 
put of the present existing simulation programs. Us- 
ing some specific modules, the user can easily write 
an interface which plugs into the 'Shower Track' the 
shower details event by event. 

The 'Shower Track' can manage all typical vari- 
ables (number of charged particles, number of elec- 
trons, energy loss..) representing the shower during 
its longitudinal development. The 'Shower Track' 
has a block of general information and a list of 
'Shower Step' objects along the EAS longitudinal 
development. 

Each 'Shower Step' object contains an informa- 
tion about the shower longitudinal development cor- 
responding to the track length dL, which may also 
vary from event to event. For each step there is also 
the possibility to handle histograms representing the 
local distributions of particles in the shower accord- 



community, to compute transmission in atmosphere. 
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ing to their energy, distance from the axis, direction 
in space, arrival time, and taking care of preserv- 
ing some of the most important correlation between 
couples of such physical quantities. 

At present it exists inside ESAF an interface to 
UNISIM, a preliminary version of interface to COR- 
SIKA, and an interface to CONEX [31,32]. 

In addition to external program interfaces, 
ESAF can also fill event by event the longitudi- 
nal description inside the 'shower object' using 
several parametrizations (Greizen [33], Greizen- 
Ilina-Linsley (GIL) [19], Gaisser-Hillas [34]) tuned 
on data and Monte - Carlo outcomes. All these 
parametrizations need the depth of the first inter- 
action point (Xi). In ESAF it can be fixed or sam- 
pled according to Nuclei- Air cross section data [35] 
or randomly. Such fast algorithms are very useful 
during the phase of detector optimization where 
the exact representation of the shower to shower 
fluctuations is less important. 

When some important information is not present 
in a given external simulation program, ESAF can 
decide either to skip or to replace using some of 
the predefined parametrizations implemented in the 
code. Several different parametrizations are avail- 
able in ESAF for energy distribution of particles in- 
side the shower [36,37,38], particles angular distri- 
bution [39], and particles lateral distribution [33,40] 
as a function of the shower age. 

Most of the results presented in this paper are 
achieved using the GIL parametrization for the lon- 
gitudinal development of the number of charged par- 
ticles N(s), as a function of the shower age s: 



3.2. Cherenkov and fluorescence light simulation 

Relativistic charged particles (mainly electrons 
and positrons) in extended air showers are at the ori- 
gin of both fluorescence and Cherenkov light emis- 
sions. Generation of fluorescence and Cherenkov 
photons has been included in ESAF. The shower 
simulation module provides for each 'Shower Step' 
the number of electrons N e (s), the electron energy 
distribution and the lateral distributions. The light 
generation module uses them to compute bunches 
of photons, as explained in section 3.2.3. This ap- 
proach is dictated by the huge number of photons 
produced which makes impossible to follow the fate 
of each individual photon. 

3.2.1. Fluorescence light simulation 

The number of photons d7Vfl uo created along dL 
by the N e (s) EAS electrons is: 



dL 
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(s) is the normalized electron energy 



spectrum, FY^(E, P, T) is the fluorescence yield (in 
photons/m) at wavelength A for an electron of en- 
ergy E. P and T are the local pressure and tem- 
perature considered at the middle of the dL track 
segment. 

The fluorescence yield dependence on the local 
temperature and density conditions is considered ac- 
cording to the empirical relation : 
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where E is the energy of the primary particle, 
A the atomic number of the primary, X the depth 
considered, X\ is the depth of the first interac- 
tion, Xq = 37.15 gem -2 (air radiation length), 
E c = 81 MeV (critical energy), a = 1.7, b = 0.76 
and Ei = 1.45 GeV. These values are chosen to 
match CORSIKA-QGSJET results [41]. 



i + b xP Vt 



where p is the atmosphere density, in kg-m -3 , and 
T is in K. The experimental values of the two co- 
efficients A\ and B\ are provided with yield mea- 
surements from Nagano measurements [43] or Kaki- 
moto ones [44]. In the Kakimoto case, the fluores- 
cence spectrum is completed with data from Bunner 
PhD thesis [45]. In ESAF, the choice to carry out 
simulations with one of these two fluorescence yield 
measurement data sets is provided by the dedicated 
configuration file. 

The yield is assumed to be proportional to the 
energy loss ^ (X is the slant depth, in g-cm -2 ) in 
the atmosphere, thus FY^ is given by: 
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FY£ (E, P, T) = FY£ (E iel , P, T) 



dE (p 
dX V^r< 



f) 



(3) 



In ESAF, the Berger-Seltzer [46] formula is used 
to compute the energy loss ^(E) directly from E. 
Fluorescence yield FY^ (E re f , P, T) at a fixed 
energy L re f is determined at T = 288 K and 
P=1013hPa either from [43] or from [44,45], 
according to the chosen configuration, with re- 
spectively L ref = 0.85 MeV and L re f = 1.4 MeV. 
Examples of computed fluorescence yield spectra 
are presented on Fig. 4. In ESAF the normalized 
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Figure 4. Fluorescence emission spectra at P = 1013 hPa 
and T = 15° C computed using Nagano measurements [43] 
(black line histogram) and Kakimoto ones [44] completed 
with Bunner rays [45] (red filled histogram). The main flu- 
orescence rays are between 330 and 400 nm, and the three 
most intense rays are at 337, 357 and 391 nm. 

electron energy spectrum is provided by the shower 
module (see section 3.1). Using equation (3), label- 
ing FY L the integrated fluorescence yield over the 
wavelength spectrum, and computing the averaged 
electron energy loss, 



<2f>. 
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we can derive from (2) the following formula imple- 
mented in ESAF: 



dN, 



fluo 



dL 
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dE \ 
dX /s 



(4) 



dx(^ref) 

The resulting total fluorescence yield (Fig. 5) 
varies weakly along the EAS development (the 
mean yield is between 4.1 and 4.7 photons/m). The 
photons longitudinal distribution along the shower 
axis is therefore similar to the electrons longitudinal 
profile. The total number of fluorescence photons 
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Figure 5. Fluorescence yield FY L (h) / FY L (0) integrated 
over the emission spectrum as a function of the altitude h 
(dashed blue line: from [43] FY L (0) = 3.81 photons/m, con- 
tinuous red line: from [44], FY L (0) = 3.26 photons/m). P 
and T profiles are obtained from the US Standard atmo- 
sphere model (see section 2.2). The mean altitude of the 
maximum of shower development is indicated for 3 shower 
zenith angles. The weak dependence of the yield with alti- 
tude is observed. 



isotropically emitted along the shower track is in 
the order of 10 15 for a primary energy of 10 20 eV. 
This number depends on the shower zenith angle, 
since inclined EAS starting at higher altitude have 
a longer path in the atmosphere and thus produce 
more photons compared to vertical EAS. 



3.2.2. Cherenkov light simulation 

The algorithm to simulate Cherenkov photon 
generation is similar to the fluorescence one. The 
number of Cherenkov photons dN c ^ ov originated by 
N e (s) electrons along a shower track segment dL is: 



dA^ ckov 



dL 
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1 dTVe 
N e dE 



{s)\CY L {E,n) 



dE 



where L t hres is the Cherenkov production energy 
threshold. The term A/"" 1 ^^ (s) describes the nor- 
malized energy spectrum of the electrons at the 
emission point. CY L (E,n) is the total Cherenkov 
yield integrated over [A m i n ,A max ]. The atmosphere 
refractive index n verifies n — 1 <C 1; it varies very 
weakly with A in the UV range of interest, thus 
only n variation with atmosphere density is taken 
into account. Therefore Cherenkov yield can be 
computed in ESAF using the formula: 
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CY L (E,n) = 
47TQ; (n — 1) 



1 



) ■ (6) 



where a is the fine structure constant, and 
^thres ~ -#=T- Typically A min = 300 nm and 

A /2(n-l) 

A max = 450 nm; these values can be configured. 
It has to be noticed that Cherenkov yield depends 
on local atmosphere parameters through the n — 1 
factor present in (6) and in the energy threshold 
expression. The normalized energy spectrum of the 
electrons at the emission point is provided by the 
shower generation module. 

Then, an effective yield (CY L ) (n) correspond- 
ing to the integral of CY L (E,n) over the electron 
energy spectrum is computed at each step of length 
dL and used to compute the number of Cherenkov 
photons produced along the step dL: 

^=N e (s) (CY L )f(n) . (7) 
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Figure 6. Cherenkov effective yield (CY 



L\ e 



as a function 



of the altitude, for 3 values of the shower age s. 

Contrary to fluorescence, Cherenkov yield de- 
creases rapidly with altitude (Fig. 6). Nevertheless, 
the number of Cherenkov photons produced all 
along the shower track (~ 3 • 10 15 at 10 20 eV) is 
of the same order as the fluorescence one. It varies 
also with shower zenith angle (Fig. 7). Cherenkov 
intensity reaches a maximum for shower zenith 
angle around 50°. In case of larger zenith angle, 
development of showers occurs at higher altitudes, 
where the Cherenkov yield is weaker. Quasi- vertical 
showers may reach Earth ground before their devel- 
opment ends, and Cherenkov photon production is 
not complete. 




70 80 

B (deg) 

Figure 7. Number of Cherenkov photons Af c k ov (#)/iV c k ov (0 o ) 
produced by an EAS (proton, E = 10 20 eV) as a function of 
the shower zenith angle 0. 

The Cherenkov emission angle C , as the yield and 
the energy threshold, is related to the atmospheric 
refractive index n and therefore correlated with al- 
titude. In the atmosphere and for shower relativist ic 
electrons, C < 1.3°. In order to obtain the angular 
distribution of the emitted Cherenkov photons, one 
has to convolute shower electron angular distribu- 
tion with C . In ESAF, the formula established by 
Elbert et al. [47] is used: 



d7V 7 
d#T 



C 



■ exp 



with #n 



0-85 £7-™° 



'7 ^0 

where C is a normalization constant such has 
Jo° ^7 = Ph°tons emitted at high alti- 

tude are beamed close to the shower axis: at 12 km, 
< 1 >~ 4°, at 4 km, < 6 1 >~ 6°. 

3.2.3. Bunches of photons 

The number of photons produced by a typical 
shower is huge, of the order of 10 15 . This fact for- 
bids to follow a direct approach, and simulate the 
fate of each individual photon. For this reason, the 
light simulation is done by introducing the concept 
of bunch, i.e. an object dedicated to EAS photon de- 
scription. 

Bunches of photons are created along the EAS 
axis: the shower longitudinal distribution is split in 
steps of length dL, set by the user in the steering 
configuration file; typical values are between 1 and 
10 g-cm -2 . At each step, two bunches are formed, 
one for the fluorescence emission, the other for the 
Cherenkov one, according to equations (4) and (7) 
respectively. Each bunch is characterized by a mean 
position value, a mean direction according to the 
computed angular distribution, and a creation time. 
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Specific distributions are associated to each bunch of 
photons, according to the emission type: wavelength 
spectrum, longitudinal distribution inside the step, 
lateral and angular distributions (w.r.t the shower 
axis). So, when a photon has to be randomly sam- 
pled from a bunch to be propagated, its character- 
istics can be retrieved from each distribution asso- 
ciated to the bunch. 

3.3. Photon propagation from EAS to the 
space-based detector 

As mentioned above, more than 10 15 UV photons 
are produced by a 10 20 eV EAS. Before reaching 
the telescope, these photons may interact with at- 
mosphere constituents or Earth ground: the physics 
processes to be considered are absorption and scat- 
tering. The tracking of such an amount of photons 
through the atmosphere is obviously impossible with 
"full" Monte-Carlo algorithms. In ESAF, two differ- 
ent approaches have been designed, each of them be- 
ing adapted to specific simulation needs. In the fol- 
lowing, the different types of interactions which can 
interfere with photon transfer towards the space- 
based detector are outlined, as well as their imple- 
mentation in ESAF. Then the two propagation al- 
gorithms are described. 

3.3.1. Scattering and absorption processes 

Photons emitted in the direction of the telescope, 
and that are neither absorbed nor scattered during 
their path arrive at the entrance pupil. To evalu- 
ate this fraction of light straightly transmitted, the 
transmission T = e~ 6 is evaluated from the interac- 
tion cross section (or from associated quantities such 
as the optical depth and the interaction length) re- 
lated to the different extinction processes, and tak- 
ing into account each photon characteristics (wave- 
length, position, ...). EAS photons that are not emit- 
ted in the direction of the detector may be scattered 
later in this direction. To quantify this effect, one 
needs to know as well the angular distributions of 
the scattering processes, represented by their phase 
function. To simulate the radiative transfer, both 
optical depths and phase functions should be com- 
puted in ESAF. 

When the scattering particle size is much smaller 
than the photon wavelength, Rayleigh theory is 
used, as for scattering on air molecules. When 
the scattering particle size is of the same order or 
larger than the photon wavelength, the scattering 



is described by Mie theory [48]. In the case of a 
space-based detector, this latter type of description 
is required to treat photon propagation in clouds 
and in aerosols. 

3.3.1.1. Rayleigh scattering In ESAF, the optical 
depth (5 ray i at a given wavelength A (in nm) is com- 
puted with the formula used in L0WTRAN7: 



*rayl 




where X is the atmospheric slant depth along the 
photon path. The A -4 dependence, at the first order, 
of Rayleigh cross-section is dominant here. The term 
in A 2 is partly due to the dependence of the cross- 
section on the air index. Studies made with ESAF 
show that vertical transmission between ground and 
top of atmosphere is around 50% at A = 337 nm, 
30% at 300 nm, and 80% at 450 nm. Simulations 
have been performed with several different density 
profiles: results show that air density profile varia- 
tions have very small impact on associated vertical 
transmission value (5T/T < 2%). 

The normalized phase function <£ ray i of Rayleigh 
scattering is given by: 

*rayl(0,p)=^(l+COS 2 0). 

Rayleigh scattering is independent of the azimuth 
angle ip and almost isotropic: the phase function 
varies by only a factor of two between 0° and 180°, 
that is very weak compared to phase functions asso- 
ciated to other scattering processes (see below, and 
Fig. 8). 

3.3.1.2. Mie scattering Below an altitude of 
around 8 km, clouds are made of tiny water droplets. 
Due to the spherical shape of droplets, the original 
Mie theory can be applied, provided knowledge of 
the distribution of the scattering droplets radius 
n(r). At a given A, n(r) determines the phase func- 
tion. In UV range, clouds do not absorb photons, 
but scatter them. Scattering can be considered as 
independent of wavelength in the range 300-450 nm 
[45]. The same behaviour is observed for high al- 
titude cirrus clouds (cirrostratus) composed of ice 
crystals. In ESAF a cloud is described as a homo- 
geneous layer, with a finite vertical height (TOVS 
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cloud). Cloud databases provide integrated vertical 
depth A vert . Laying from 1.2 to 3.3, the mean value 
of A ver t is around 2. For A ver t — 2, the transmission 
of a vertical UV light beam crossing the cloud layer 
is - 15%. 

Phase function values for low altitude clouds, ob- 
tained from one particular cumulus model [49], are 
stored in tables for A = 450 nm, as well as those 
for cirrostratus, for A = 500 nm [50]. Contrary 
to Rayleigh scattering, Mie scattering in clouds is 
highly anisotropic, and very forward peaked. Phase 
functions are compared in Fig. 8. 
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Figure 8. Cloud phase functions (cumulus and cirro- 
stratus) compared to Rayleigh one. P(9) is equal to 
27r$(6>,(^) = 2>k<$>{6). 



Concerning aerosols, the three main models (ru- 
ral, urban, maritime) provided by L0WTRAN7 have 
been introduced in ESAF. These models provide 
the dependence of the extinction coefficient with 
altitude and wavelength, specifying the respective 
part of scattering and absorption processes. These 
aerosol layers are located between and 2 km alti- 
tude. More details can be found in [23] (in french) 
or [51]. The aerosol layer opacity is specified via the 
visibility. The smaller is the visibility, the greater is 
the extinction coefficient. Default visibility values in 
L0WTRAN are 5 km and 23 km, corresponding to a 
vertical transmission of 75% and 20% respectively. 
To describe scattering angular distribution, ESAF 
is supplied with phase functions from L0WTRAN7 for 
the three mentioned models. 



3.3.1.3. Ozone absorption In the UV range of in- 
terest, the molecular absorption in atmosphere is 
dominated by ozone (O3) absorption. At a fixed al- 
titude h, the ozone absorption length L a bs,o 3 is ex- 
pressed by 

— — — = Po 3 (h) cr abS5 o 3 (A, h). 

^abs,0 3 

Ozone absorption cross section (cr a b s ,o 3 ) data pro- 
vided by L0WTRAN have been introduced in ESAF. 
The dependence of cross-section on local tempera- 
ture is taken into account, entailing a dependence on 
altitude. po 3 (h) is determined from available ozone 
profiles, and air density. Ozone exists mainly be- 
tween 10 and 30 km, with a maximum absorption 
around 20 km. It absorbs the UV light below 330 nm, 
therefore it has a weak impact on fluorescence trans- 
mission, since the three main fluorescence lines are 
above the wavelengths affected by ozone; it affects 
mainly Cherenkov light. 

3.3.1.4. Ground reflection The reflection on the 
Earth surface is a crucial point for Cherenkov sig- 
nal detection. The Earth albedo, in UV range, has 
been measured by ULTRA [52] experiment and is 
less than 10% [53]. Moreover, reflected photon an- 
gular distribution depends on the reflection type, 
which can be Lambertian (reflected light is scat- 
tered isotropically by a Lambertian surface) or spec- 
ular (highly anisotropic, as on the ocean surface). 
Both types have been implemented in ESAF. To 
describe specular reflection, a simple model, repro- 
ducing qualitatively air-born measurements of the 
ocean BRDF (Bidirectional Reflectance Distribution 
Function) [54] with dependence on wind speed, has 
been introduced. For a wind speed of 9 m/s, specu- 
lar peak amplitude is 5 times lower than for a wind 
speed of 2 m/s. If specular reflection is dominant, re- 
flected light intensity is not isotropic in # r , and can 
be zero on a large zenith angle range (see Fig. 9). 

3.3.2. Propagation algorithms 

Below are described two original propagation al- 
gorithms designed and implemented in ESAF. 

3.3.2.1. "Bunch" algorithm The algorithm has 
been designed to handle mainly three components 
of the light signal received by the telescope. The 
first component, called 'Direct Fluorescence', is the 
fluorescence light reaching directly the lens without 
undergoing any interaction in the atmosphere. The 
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Figure 9. Simple model for light reflection on ocean. The 
zenith angle #i of the incident light beam is 20°. The re- 
flectance is computed in the plane containing the incident 
direction, as a function of the reflected light zenith angle 6 r , 
in case of both lambertian or specular reflection (for 2 wind 
speed values u). 

second component, called 'Reflected Cherenkov', is 
the ground reflected Cherenkov light, having inter- 
acted only once along its path from emission point 
to the telescope. The third component, called 'Air 
Scattered Cherenkov', represents Cherenkov pho- 
tons which have been scattered only once, on air 
molecules. 

The extreme low value of the detector solid an- 
gle (Q ~ 10 -11 sr for EUSO telescope shipped on 
ISS) is of major importance: when considering pho- 
tons emitted or scattered within the detector solid 
angle, 'single' photons can be sampled from the re- 
lated bunches, and be transmitted one by one to the 
telescope. 

In the case of fluorescence photons directly emit- 
ted towards the detector, single photons are sam- 
pled from bunches created at fluorescence emission 
point. Then, they are distributed in time, position 
and wavelength, according to distributions associ- 
ated to the bunch they originate. For each photon, 
the probability that it is transmitted up to the tele- 
scope is computed taking into account all possible 
interactions in atmosphere. 

Concerning the two other components, they are 
both simulated by the following algorithm. Since 
Cherenkov radiation is beamed in the forward di- 
rection, it is possible to propagate each Cherenkov 
photon bunch step by step, until it reaches ground 
or exits atmosphere. At each step, two operations 
are performed: 

(i) The number of Cherenkov photons scattered 



once on air molecules, and directed toward the 
telescope after the interaction, is computed. 
A corresponding number of single photons are 
therefore created. Their final transmission to 
the telescope from their scattering position is 
then evaluated : if they undergo another inter- 
action along their path to the detector, they 
are lost. 

(ii) The wavelength spectrum of the propagated 
bunch is convoluted by transmission values 
computed along the step. Then, the possibility 
to pursue the bunch propagation is evaluated. 
When a Cherenkov bunch reaches the ground, the 
number of photons reflected toward the detector is 
computed. A corresponding number of single pho- 
tons are created and their transmission to the tele- 
scope is evaluated : if they undergo another interac- 
tion, they are lost. 

In presence of aerosols or clouds, the bunch algo- 
rithm takes into account their effects on the trans- 
mission values, considering the photons scattered on 
these media as lost. 

The bunch algorithm is useful to perform fast sim- 
ulations (see Fig. 10). But such an algorithm can 
not handle properly single scattering of fluorescence 
photons, and multiple scattering of both fluores- 
cence and Cherenkov photons. A mandatory point 
in EAS analysis is to evaluate the detailed contri- 
butions of scattered photons to the detected signal, 
as this component may introduce systematic uncer- 
tainties if not taken into account during reconstruc- 
tion procedure. To achieve this, a more sophisticated 
algorithm has been designed to deal with single and 
multiple scattering of the light produced by EAS. 

3.3.2.2. "Reduced" Monte- Carlo algorithm First 
of all, the simulation of 'Direct Fluorescence' com- 
ponent does not raise any problem, and it is treated 
in the same way as presented in the previous sec- 
tion. Then, the aim is to simulate fluorescence 
and Cherenkov light propagation by performing 
a Monte-Carlo ray-tracing algorithm, propagating 
photons one by one. Let us consider photons that 
reach the telescope after exactly n scattering pro- 
cesses. To simulate the transfer of each of them, an 
ideal algorithm has to follow the steps described 
below, as illustrated on Fig. 11. 

(i) Photon position, direction and wavelength 
have to be randomly generated, taking into 
account related distributions defined during 
light production process (see section 3.2). 
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Figure 10. Examples of time photon distributions, for a 
shower with E = 10 20 eV, = 60°, and a footprint at the 
telescope nadir - The different contributions are superim- 
posed: fluorescence photons transmitted directly to the tele- 
scope, Cherenkov photons scattered once by the atmosphere 
molecules, and Cherenkov photons reflected by the ground 
in case of clear sky (top panel) or scattered by the cloud in 
case of cloudy sky (bottom panel). 

(ii) Interaction type and position are randomly de- 
termined according to the different physical 
processes able to alter photon path. If the pho- 
ton is absorbed or if is going out of the atmo- 
spheric medium before interacting, it is con- 
sidered as lost. 

(iii) Previous step is repeated on the scattered pho- 
ton, until its loss or until n-th scattering. 

(iv) At the nth scattering, a random sampling is 
performed considering the probability : 



p = <s> n {o^)n det 



(8) 



(v) 



in order to determine whether the photon is 
scattered toward the telescope or not. <£ n (0, tp) 
is the phase function of the n-th scattering, 
and fldet the solid angle associated to telescope 
pupil entrance. If the photon is not scattered 
toward the telescope, it is rejected. 
Then, transmission from the n-th scattering 



position to the detector is computed. If the 
photon undergoes another interaction it is lost; 
if not, it reaches the telescope lens. 

"full" MC: N photons 
"reduced" MC: N re d photons 
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Figure 11. Scheme of the Monte-Carlo (MC) algorithm prin- 
ciple to simulate photons scattered exactly n times. 

As previously stated, with the ideal algorithm de- 
scribed above, it will be unaffordable CPU time con- 
suming to propagate each of the ~ 10 15 UV pho- 
tons produced by an EAS initiated by a primary of 
10 20 eV (see section 3.2). About 10 4 photons reach 
the telescope. The high decrease between these two 
numbers is due to the solid angle. If this reducing fac- 
tor could be applied before photon tracking, a ray- 
tracing simulation should be possible. The key point 
comes from the particular situation of a space-based 
detector : it stands outside the scattering medium. 

Above an altitude of 30 km, atmospheric density is 
very weak, and photon scattering can be neglected. 
As a consequence, one can define a minimum dis- 
tance D m i n below which photons cannot be scat- 
tered. This distance and the associated solid angle 
l] max = ^det/^min are snown on Fig. 12. In case of 
a detector shipped on ISS, its altitude is ~ 430 km, 
thus f^max = 3 • 10 _11 sr. 

The probability P defined by the equation (8) can 
be written as the product of two independent factors 
P = Pi x P 2 defined as: 



Pi 



$ n (0,(p)Sl det 
^max^max 



max^ L max 



(9) 
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Figure 12. Minimum distance D m i n between the top of at- 
mosphere (considered here at an altitude of 30 km) and the 
telescope and the associated largest solid angle Slmax- 

$ max is the maximum over all phase functions in- 
volved in simulation. By normalising the phase func- 
tion 3> n (0, (/)) of the nth scattering (which type is a 
priori not known) by $ max , the condition Pi < 1 is 
insured. £l max is small enough to insure P2 < 1. Ran- 
dom draw of P probability is split into two indepen- 
dent draws of P\ and P2. Since P2 has been set up 
to be independent of all other stochastic processes 
involved in photon propagation, random computa- 
tion of P2 can be performed at the very beginning of 
simulation chain. Moreover, considering that it can 
be identically applied to all photons, only one evalu- 
ation is performed, according to a binomial law with 
(TVo, P2) parameters. No is the total number of pho- 
tons produced by the shower, Nq ^> 1 ^> P 2 . Then 
the reduced number of photons to propagate iV" re d 
when simulating the nth scattering, is given by: 

7V red = PoisSOn(7Vo$max^max) , (10) 

where the Poisson distribution is the binomial dis- 
tribution limit, in case of a large sample size and 
low probability. The described algorithm is then ap- 



plied, taking probability P\ instead of P in the 4th 
step. All interaction types described in section 3.3.1 
- molecular absorption, scattering on air molecules, 
aerosols, clouds, and ground - are treated by this al- 
gorithm. 

The treatment presented above is valid for a fixed 
scattering order. To obtain the contributions of dif- 
ferent scattering orders n and p, the described al- 
gorithm has to be performed for one set of gener- 
ated photons N^J , until the nth scattering, then for 
another set until the pth scattering, in a full 

independent way. Before computation, a maximum 
scattering order n max has to be fixed. The contribu- 
tion of the scattering processes to the signal is the 
sum of the photons scattered exactly 1 time, 2 times, 
3 times, n max times. "Zero" order corresponds to 
fluorescence photons produced in the detector solid 
angle, and transmitted without interacting up to 
the telescope. Considering the first scattering order 
in clear sky conditions, the different contributions 
come from: 

• zero order 

- fluorescence photons straightly transmitted 

• first order 

- Cherenkov photons reflected on ground 

- Cherenkov photons undergoing molecular scat- 
tering 

- fluorescence photons reflected on ground 

- fluorescence photons undergoing molecular 
scattering. 

The main advantage of this algorithm is to be 
able to simulate radiative transfer of photons emit- 
ted by any type of light source, taking into account 
the multi-scattered photons at any order; but its 
main limitation is the computing time, which can 
be high, despite the applied reducing factor. The 
computing time is strongly dependent on maximum 
scattering order n max , as well as on number of pho- 
tons to propagate iV r ed- The latter is related to the 
maximum solid angle, fixed by detector setup, to 
the number of emitted photons, related to shower 
parameters (22, 6) and to the maximum value of the 
phase function. $ max can vary from 0.3 in case of 
clear sky condition to 1000 in case of cloudy condi- 
tion. Simulation runs show that computing time is 
not affordable in the latter case. 

To bypass this limitation, the cloud phase func- 
tion involved in the last scattering process is ap- 
proximated. Let us highlight that the cloud phase 
function involved in previous scattering processes is 
not approximated ; it concerns only the last scatter- 
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ing process, if it occurs within a cloud. Two kinds of 
approximation have been studied. The first one uses 
an analytical form, the Double Henyey-Greenstein 
function [55] . The second one uses an averaged phase 
function peak over a limited range of 10°. 

In the current state of ESAF development, aerosol 
phase functions are treated in a simpler way, us- 
ing the analytical approximation for every scatter- 
ing process occurring on aerosol layer. 
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Figure 13. LIDAR shot at 1064 nm from ground to 1 km 
distant cumulus cloud, 300 m deep and with an extinction 
coefficient of 17.25 km -1 (no atmosphere contribution, cloud 
only). Contribution of the multiple scattering component 
for a receiver field of view of 10 mrad (top) and 1 mrad 
(bottom). ESAF results are obtained in simulating scattering 
process until order 8. In both cases the ESAF results are 
superimposed to MUSCLE ones, showing that they are in 
very good agreement. 

3.3.3. Validation tests 

To proof that both propagation algorithms are ef- 
ficient to properly simulate radiative transfer, vali- 
dation tests were carried out. First, the propagation 
in clear sky conditions of photons emitted by a given 
shower (fixed energy and zenith angles) has been 
simulated by the two designed algorithms. Results 
concerning Cherenkov photons propagation (molec- 
ular single scattering, ground reflection) were com- 



pared. The time distribution of photons at telescope 
entrance pupil and their last scattering position in 
FoV are similar, clearly proving the consistency of 
the two algorithms. 

In a second step, studies have been undertaken 
to test the validity of the most efficient algorithm, 
namely the "reduced" Monte-Carlo algorithm. The 
aim is to compare it to simulation tools developed 
within the atmosphere physicists community. One 
study, handled by MacKee et al. [56], concerns 
cloud albedo change with optical depth, predicted 
by a Monte- Carlo simulation and checked experi- 
mentally. Results got with ESAF are in agreement 
with those obtained by MacKee et al. The sec- 
ond set of tests were done to reproduce studies of 
the MUSCLE (Multiple Scattering Lidar Exper- 
iments) working group [57]. The aim is to predict 
LIDAR shot response after multiple scattering in 
clouds. The results obtained with ESAF are com- 
pared to the MUSCLE ones: a good agreement 
has been found (see Fig. 13 and reference [23] for 
more details) that actually proves the validity of 
the "reduced" Monte-Carlo propagation algorithm 
implemented in ESAF. 

Let us specify that the two kinds of cloud phase 
function approximations have been compared dur- 
ing these tests. It has been shown the use of an av- 
eraged peak is the best solution, and will be used 
further in this paper. 

With this propagation package in ESAF, gener- 
ated fluorescence and Cherenkov lights are trans- 
fered up to the space-based detector. Thanks to the 
designed algorithms, careful studies of the signal at 
detector pupil entrance are possible dealing with 
different atmospheric situations: clear sky, aerosol 
layer, clouds. 

3.4. Detector response simulation 

The ESAF detector simulation is the part of the 
framework which is more tightly connected to the 
EUSO design. Most of the implementation is based 
on a monocular cylindrical telescope with a photo- 
multiplier (PMT) based photodetector. The main 
scheme of the detector is shown in Fig. 14. Never- 
theless, the ESAF detector simulation structure re- 
tains large degree of flexibility and can be adapted 
to a different detector concept. 

The telescope simulation part is divided in 2 main 
modules: 

(i) The Monte-Carlo photon tracking engine 
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through the optical system, 
(ii) The electronics and trigger response. 

3.4.1. Optics and focal surface 

The detector simulation starts from the photons 
in front of the entrance pupil. One by one the pho- 
tons are propagated through the optical elements of 
the detector; in Fig. 14 the baffle, the optical system 
and the focal surface are visible. Three simulated 
optical systems are provided within ESAF: 

(i) The EUSO Phase A baseline optics, a double 
Fresnel system. The tracking through this op- 
tics is simulated by tracking each individual 
photon entering the pupil up to the focal sur- 
face. A detailed description of the fine struc- 
ture of the Fresnel lenses allows to reproduce 
and study the optics distortions. As shown in 
Fig. 15, a large fraction of the light in the Point 
Spread Function (PSF) of the double- Fresnel 
optics is spread all over the focal surface: the 
so-called stray-light. A commercial program 
(CodeV, [42]) was used to simulate outside 
ESAF the exact ray tracing properties of the 
lenses and the result of this CodeV simulation 
was used to create an accurate parametriza- 
tion of the PSF, of the amount of stray light 
and of the aberrations. Within ESAF, this 
parametrization is used to trace the photons. 
We underline that the light that contributes 
to the signal and, particularly, to the trigger, 
is the one that is contained in a circle of the 
size of the focal surface pixel. Fig. 16 shows 
that in a typical pixel of about 5 mm diame- 
ter (EUSO design), about 60% of the photons 
contribute to the trigger, while the other 40% 
is lost as stray light. Fig. 16 shows also the 
(weak) dependence of this number as a func- 
tion of the angle that describes the position of 
the shower in the field of view (see Fig. 23). 

(ii) A second multi-purpose optical system, whose 
characteristics are implemented through a 
parametrization of the optical response (PSF, 
stray light, aberrations) and are defined by 
a set of configuration files. Highly flexible, 
it allows an effective fast simulation without 
tracking individual photons. It is particularly 
useful for studies where the performances of 
the optics have to be tuned according to the 
overall requirements, or as a reference to com- 
pare the performances of the full simulation, 
or to speed up the simulation for large statis- 



tics studies. The parametrization and the con- 
figuration files were obtained by running the 
full simulation. Most of the work in this paper 
was obtained with this faster implementation 
but we have checked that the results obtained 
with the full simulation are always consistent, 
(iii) A Geant4 based implementation of the optics 
and of the ray tracing (see next paragraph). 
This Geant4 implementation is new and was 
not used in this work. 
Ideally the efficiency of the optics scales with the 
angle in the FoV, 7, as cos 7; for the EUSO optics 
the mean efficiency is of about 45%. To this number, 
another 60% light loss must be added because of the 
optical aberrations and stray light, as determined 
from the simulations. 



Focal surface Optical System 




Figure 14. Main scheme of the detector. 



The focal surface module implements the geome- 
try of the EUSO photodetector. The two available 
layouts of the focal surface are illustrated in Fig. 17. 
The focal surface is based on the modular struc- 
ture of small autonomous units called elementary 
cells [58], each one being a matrix 2 x 2 of multi- 
anode photomultipliers. An optional optical adapter 
can be placed in front of the PMTs (efficiency of 
about 95%). The focal surface elementary cell lay- 
out geometry is loaded at runtime from the layout 
configuration file. The simulations show that the fo- 
cal surface, in the polar layout, has a "filing factor" 
efficiency of ~ 90%. The tracking engine propagates 
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Figure 15. An example of the irregular PSF of the Fresnel 
optics generated using an incident flux at an angle 7 = 25 deg 
in the FoV. (a) Distribution of photons on the whole focal 
surface, (b) Zoom on the spot of the previous picture. The 
red circle corresponds to a bucket of 5 mm of radius. 



the photons up to the PMT photocathodes and con- 
verts them into an electronic signal, stored for the 
electronics simulation. 

Any PMT (or optical adapter) that fulfills the 
elementary cell specifications and constraints can 
be easily added to the focal surface simulation. 
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Figure 16. Fraction of energy (fraction of photons) contained 
in a circle as a function of the circle diameter. 




(a) 



(b) 



Figure 17. Focal surface layouts in ESAF. (a) Cartesian 
layout, (b) Polar layout. 

Currently two models 6 are available, the 36-pixel 
Hamamatsu R8900 and the 64-pixel Hamamatsu 
R7600. 

In the Fig. 18a and 18b there is an example of how 
an EAS appears on the telescope focal surface. 

3.4.2. GEANT4 simulation of optics in ESAF 

A detector simulation chain with the possibility 
to trace particles through the detector using the 
GEANT4 package [59] has been implemented in 
ESAF. We have kept the original way of ray-tracing 
for checks and compatibility. This upgrade makes 
ESAF much flexible for a simulation of any detector 
and not necessarily an EUSO-like detector. As two 
examples we implemented an EUSO like detector 
with two Fresnel lenses and a diflxactive surface 
between them and a TUS-like detector [60] with a 
Fresnel mirror. Both examples are shown in Fig. 19. 



6 In the PMT models ESAF takes into account the quantum 
efficiency as a function of the wavelength, the collection 
efficiency, the gain, the cross-talk between pixels, the PMT 
geometry, the time width of the signal and the dark noise 
rate. The typical PMT efficiency is 25% (quantum) plus 70% 
(collection). 
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focal surface. No background is simulated, (a) 2D view. The 

small squares are the pixels, the big ones the PMT. (b) 3D Figure 19 EUSO Fresnel lens (top) an TUS Fresnel mirror 

view including also the time counted in Gate Time Units (bottom) detectors built with Geant4. Fluxes of photons 

(GTU) of 2.5 fjs. The x and y axis are the xy coordinates of focused on their respective focal planes are also displayed, 

the pixels while the z axis is the time from the first signal 
photon counted (in GTU). On the xy-plane is plotted the 
profile of the pixels with signal. 



Almost all segments of EUSO Fresnel lenses have 
parabolic surfaces. Every segment is done with 
boolean operations on the basic shapes (G4Tubs 
and G4Paraboloid). More complex surface shapes 
were approximated by cones. Approximation pre- 
cision can be chosen on demand. For the current 
simulation it is equal 0.1 jiva. A piece of a EUSO 
Fresnel lens implemented with GEANT4 is shown 
in Fig. 20. 

The following standard optical processes are de- 
fined in the GEANT4 implementation: Fresnel re- 
fraction, Fresnel reflection, total internal reflection 
and absorption. In addition to the standard optical 
processes, ESAF uses a "special" process parametri- 
sation: the diffraction on the second Fresnel lens. In 




Figure 20. A piece of a EUSO Fresnel lens simulated with 
help of GEANT4 
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fact, the second lens of the EUSO detector has two 
surfaces, one with Fresnel structure and the second 
one with a normal diffractive surface. The role of the 
diffractive surface is to reduce the chromatic aber- 
rations. For the simulation of the diffraction pro- 
cess on the second lens we use a parametrization. 
The momentum direction correction depends on the 
material refraction index and on transverse compo- 
nent of position on the lens surface. To implement 
it in the Geant4 simulation we have created a spe- 
cial class that inherits from G4VDiscreteProcess. 
We have added a thin vacuum layer just after the 
flat lens surface, and have assigned proper material 
parameters. Thus, rays going out of the lens are af- 
fected by ordinary boundary processes, and, if not 
reflected back, they go through that thin vacuum 
layer and are therefore handled by the new class. 

The GEANT4 implementation of the simulation 
was not used in the work described in this paper. 
For this reason, we do not give a detailed account 
of GEANT4 - ESAF capabilities. However this ca- 
pability will be very important for future detector 
simulations and it will be documented in a future 
work. 



3.4.3. Electronics and trigger 

The electronics simulation converts all the PMT 
signals into front-end electronics (FEE) counts, sim- 
ulates the front-end chip response and applies the 
trigger patterns. 

In the front-end the PMT signals are sorted by 
time, divided in time units called Gate Time Units 
(GTUs) and counted. In the baseline configuration 
of EUSO the duration of a GTU is 2.5 /is. An am- 
plitude threshold and the double hit resolution of 
the front-end chip are taken into account as well, 
according to the specifications of the EUSO Phase 
A design report [16]. The electronics has a typical 
efficiency of about 90%. 

The atmospheric background contribution is 
added directly after the front-end simulation as 
additional FEE counts. The motivations for this 
solution are described in detail in section 3.5. 

The trigger in ESAF is implemented as a set of 
'Trigger Engine' objects, one for every trigger algo- 
rithm. Several engines can be applied to the event, 
giving the possibility to compare their efficiency on 
an event-by-event basis. The specific trigger algo- 
rithm used for the studies presented in this paper is 
described in the section 3.4.4. 



3.4.4. The level- 1 trigger algorithm: Tracking 

The image of a developing EAS on the focal sur- 
face is a bright spot moving on a straight line whose 
brightness rises up to the maximum and gradually 
fades away. At a given time, the spot is mainly con- 
tained in a single pixel, as required by the optics de- 
sign. Integrated over the development, the shower 
appears as a track in space and time whose length 
and brightness depend on the primary energy and 
incident angle but also on the EAS position in the 
FoV. 

In a space-based experiment, only those EAS that 
produce a track on the focal surface, standing out 
from the surrounding nightglow background, can be 
detected and provide any useful information. With- 
out the track structure visible, the EAS signal is not 
distinguishable from the background fluctuations. 

Therefore two conditions can be established for 
an EAS to be detectable: 

Brightness: The signal in a given pixel and GTU 
must be well above the average background 
counts. 

Contiguity: The bright, or active, pixels must be 

contiguous in space and time. 

At pixel level, to set a "brightness" requirement, 
a background dependent threshold M can be estab- 
lished: 

M{k) = n h + k^ , (11) 

where is the mean number of background counts 
per pixel per GTU, and k a coefficient depending on 
the requested background rejection efficiency. Since 
the random background is a poissonian distribution 
with mean nt>, k can be regarded as a number of 
standard deviations. This parametrized threshold is 
intrinsically adaptable: for k fixed, can be con- 
tinuously measured and the threshold can be auto- 
matically updated. 

The contiguity criteria requests that the signal 
must form a track in space and time. Starting from 
the first active pixel, the next active one in the next 
GTU must be the same pixel or one of the nearest 
neighbors. The minimum number of consecutive ac- 
tive pixels, the track length L, can be assumed as 
the second requirement a signal must satisfy to be 
recognized as a possible EAS candidate. 

We have studied the trigger efficiency using the 
contiguity tracking trigger (CTT) developed dur- 
ing the EUSO Phase A. It works at the Elemen- 
tary cell (144 pixels) front-end level, that is an en- 
semble of 2 x 2 closely packed multianode PMT (36 
channels each). The CTT is a simple algorithm that 
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checks the signal contiguity in space and time, as 
explained before. Thought to be integrated in the 
front-end chip as first level trigger, it can be imple- 
mented using purely combinatorial circuits with low 
cost in term of power and mass budget. Due to the 
large background non-uniformities, M cannot have 
the same value for every channel on the focal sur- 
face. Instead, we used k as one of the variables of the 
trigger. In the following we shall denote a particular 
CTT configuration with the notation [L; k]. 

The rate of fake triggers (over the whole detec- 
tor) due to the random background is easy to com- 
pute [21,22]: 

N cell P(L;n h ;X(k)) 



z/ftr [L\ k] 



^GTU 



(12) 



where N ce \\ is the number of elementary cells (1386 
in the EUSO Phase A configuration), £gtu is the 
duration of a gate time unit and 



P(L;n h ;AT(k))^9^ 



M(k)- 

£ 

i=0 



is the poissonian probability to have a random track 
of length L with a pixel threshold Af(k). 

The maximum acceptable overall fake trigger rate 
is constrained mainly by the available telemetry re- 
sources. Thus, the maximum fake CTT rate is then 
limited by the time the second level trigger further 
processes the data. The final number depends of 
course on the details on the second level trigger, but 
a realistic estimation, based on the current technol- 
ogy, sets this limit to ~ 100 Hz. 

3.5. Simulation of the random background 

The main source of background to the EAS signal 
is the atmospheric nightglow (see [61] and references 
therein). Other sources like moon light, lightenings 
and man-made lights are not constant in time. Cur- 
rently, ES AF is able to simulate only the effect of an 
uniform background. This is not a major issue since 
the main interest is to study the detector capabili- 
ties with various levels of average background. The 
mean value of the background radiance used in this 
work is B « 0.5^2 -10 12 photons m~ 2 s" 1 sr" 1 , as 
measured by different experiments (e.g. BABY [62], 
TATIANA [63]). 

Assumed uniform and isotropic on the EAS devel- 
opment timescale, it floods the entrance pupil with 
a constant flux of photons. The resulting flux on the 



focal surface is not uniform, because of the optics ef- 
fects. The average background counts per pixel per 
GTU, ^ 1X ' vai> i es from the focal surface center to 
the border with a profile that depends on the optics 
used. ^£ lx sums up the contributions of photons from 
all incoming directions and therefore the stray-light 
(see Fig. 15) cannot be neglected. 



z/£ lx (x) = p h (x) s PB d 2 



(13) 



where x is the pixel position, pb the number of back- 
ground photons per surface unit, £pd the photo- 
detection efficiency and d the pixel side. The distri- 
bution pb (x) is obtained by simulating several mil- 
lion photons with an uniform incident direction dis- 
tribution, according to the background radiance. For 
such a purpose ESAF has a special running mode 
where only the optics is loaded, optimizing the mem- 
ory and CPU consumption. 

Once this has been done, and accordingly with the 
user choice, ESAF can simulate the random back- 
ground over the whole focal surface only in those ele- 
mentary cells that have at least one pixel with signal. 
In each pixel the simulation adds a random num- 
ber of background counts for every GTU, following 
the Poisson distribution with mean = ^£ lx tGTU> 
where £qtu 18 the GTU time length. 

4. Shower reconstruction 

After the observation of one event, the following 
questions should be addressed: what are the direc- 
tion, the energy and the nature (nucleus, gamma, 
neutrino, . . .) of the cosmic ray at the origin of the 
detected EAS ? 

ESAF has a reconstruction package which con- 
sists of special modules aiming to recover these in- 
formations from an event. Let us consider below the 
underlying principles of the reconstruction software 
which consists essentially of the following modules: 
Pattern recognition: the goal of this module is to 
find the "shower track", i.e. a set of spatially cor- 
related pixels on the focal surface that describe 
a moving point as a function of time (the projec- 
tion of the EAS development on the focal surface). 
The main problem here is to disentangle the pix- 
els that contain some signal to those that are pure 
background. 

3D Direction reconstruction: the goal of this 
module is to find the direction of the primary 
particle in space. The 3-dimensional reconstruc- 
tion is possible by using both the direction of 
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the projected track on the focal surface and the 
arrival time of the photons. 
X max and Energy reconstruction: the goal of 
this module is to find the altitude in the atmo- 
sphere where the EAS reaches its maximum of 
development and the corresponding atmospheric 
slant depth, and to reconstruct the energy of the 
primary particle. The fluorescence light emitted 
by an EAS is proportional to the deposited en- 
ergy. However, the collected light is not because 
there are important corrections that depend on 
direction, position on the focal surface, and de- 
tector efficiency The X max and its variation with 
the energy are very important for primary par- 
ticle identification. Since the energy and X max 
are correlated, it is not possible to infer a precise 
determination of each variable in an independent 
way. Therefore, this module reconstructs the best 
estimate of the energy and of the X max . 

4.1. Pattern recognition 

Once a shower track has been triggered by the 
apparatus, the pattern recognition module selects, 
within every considered time interval (GTU), the 
pixels containing signal from the pixels containing 
background only. Since the background light is ex- 
pected to be randomly distributed both in direc- 
tion and time, the space-time correlation in the focal 
plane is the key tool to select clusters of pixels which 
belong to the shower track on the focal surface. 

The pixels selection is implemented in ESAF by 
a simple multi-step analysis which, within the fo- 
cal plane, searches for clusters of pixels in space and 
time separately. A following module also analyzes 
the X versus time and Y versus time correlations in- 
side the focal plane by means of a linear fit, followed 
by the exclusion of the points far away from the fit- 
ted line. A more sophisticated method using similar 
concepts, but based on the Hough transform tech- 
nique [64] is also successfully implemented in ESAF. 

4.2. Direction reconstruction 

Let us introduce a standard reference system 
(right-handed) with the origin in the center of the 
detector focal surface and the z-axis pointing to- 
wards the nadir of the instrument (see Fig. 23). 7 



7 We are assuming that there is no tilt between the nadir 
and the orientation of the instrument optical axis. 



The EAS axis is described in this reference system 
by an unit vector ft in terms of two polar angles 
and cp: 

ft = (sin 6 cos 0, sin sin (/>, cos 0) . ( 14) 

In an analogous way we define the unit vector of the 
direction that is seen by the i-th. front-end pixel in 
the FoV: 

ni = (sin 0\ ow cos (p} ov , sin 6\ ow sin $ ov , cos 0\ ow ) . 

(15) 

Therefore, for each selected pixel on the focal 
plane there is an associated direction in space 
where the fluorescence light activating such pixel is 
expected to come from. 8 

Such directions in space are expected to belong, 
within instrument resolution and neglecting the 
shower lateral extension, to the plane containing 
the shower axis and the detector itself. Using the 
information of the light arrival time in each pixel, it 
is then possible to look for the unique track in space 
associated to the measured direction-time pattern. 
Such pattern depends strongly on the direction of 
the shower axis in the detector- shower plane, but 
very weakly on the distance of the shower maximum 
to the detector because of non-proximity effect. This 
is one of the advantages of space-based observation: 
the altitude of the shower maximum may vary up to 
20 km, depending on the shower inclination, but it 
is always a small distance compared to the satellite 
typical altitude of 400 km, which means that the 
collected light changes at most by 10%. 

With a monocular detector, it is not possible 
to determine by means of a fitting procedure the 
detector-shower distance, unless the resolution is so 
good to allow a measurement of the spatial width 
of the focal surface track. The shower detector dis- 
tance can be simply fixed to the expected average 
value for each zenith angle of the incoming shower. 
However, such an information can be inferred by 
the arrival time of the Cherenkov light reflected 
on the ground, if it is available, since it is strongly 
correlated to the distance of the shower from Earth 
and detector. With the detector-shower distance 
inferred by such procedure, the fit of the direction- 
time pattern yields a better determination of the 
shower-axis direction inside the shower-detector 
plane. In case the Cherenkov light is not detected, 



8 For any detector configuration, a map that associates each 
pixel with a unique direction in the FoV is built simulating 
an uniform photon flux in the FoV through the detector. 
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the detector-shower distance can be inferred by the 
spatial width of the fluorescence light profile (see 
section 4.3). In both cases, the determination of the 
detector-shower distances and of the zenith angle is 
made more precise using an iterative reconstruction 
process. 

The first step that leads to the determination 
of the EAS direction is the reconstruction of the 
Shower- Detector Plane (SDP), i.e. the ideal plane 
that contains both the EAS track and the detec- 
tor. The SDP can be reconstructed by ESAF fitting 
the "moving" EAS image on the planes x-time and 
y-time. In this way one obtains a projection of the 
shower track on the plane tangent to the Earth sur- 
face at the telescope nadir point. The SDP is then 
the plane that contains this projection and the de- 
tector itself. 

The second step is the fit of the track on the SDP 
in order to infer the direction of the shower axis. 
In ESAF two different classes of track fitting are 
implemented [65]. 

Numerical Exact (NE) For each signal FEE 
count we compute its expected arrival time, 
^(rii,fi,a), as a function of the EAS axis di- 
rection (unknown), the pixel direction in the 
FoV (see above) and other free parameters 9 , 
a. These times can be exactly calculated using 
only geometrical and physical relations. Since we 
know the experimental arrival times, we can 
minimize the % 2 function 

A 2 

2 /o \ _ (t* -ti(ni,«,tt)) 



X 



where the o~i are the uncertainties on the mea- 
sured arrival times, via numerical minimization 
methods. From the minimization results we know 
the EAS direction. 
Analytic Approximate (AA) In this case we 
use approximate formulas, e.g. assuming that 
the EAS track has a constant angular speed on 
the focal surface, that can be analytically solved 
using linear fit procedures. 

4.3. X max and Energy reconstruction 

The principle of the method described in this sec- 
tion is to assume a functional form Fi ong for the lon- 



gitudinal development of the shower, to calculate the 
expected number of photoelectrons at the telescope 
as a function of time and to vary the parameters of 
the shower function Fi ong until an agreement with 
the measured signal is obtained. From these param- 
eters, one can deduce the energy and X max of the 
EAS. 

To calculate the number of photons emitted along 
the track and the fraction of them that reach the 
telescope, one has to know the shower track posi- 
tion in space and the local atmospheric parameters. 
Moreover to obtain the expected signal, one has to 
estimate the detector efficiency. 

In what follows 6f ov and 0f ov as functions of time 
are needed to fix the direction in space where the 
fluorescence light is expected to come from as a func- 
tion of time. These functions are obtained by fitting 
#f ov and </>f ov of each active pixel with presumably 
a signal by a linear function of time. An example of 
such fit is shown in Fig. 21. 




9 These free parameters depend on the particular track fit- 
ting algorithm. They can be, for example, the position in 
FoV, the time and the distance from the detector of the EAS 
maximum. 



GTU (us) 

Figure 21. An example: vs time (upper plot) and cf) vs time 
(bottom plot) (each GTU is 2.5 fis). 



4.3.1. Background and signal estimations 

A cluster of pixels belonging to the shower track 
is taken from the pattern recognition module to as- 
sign for each time stamp a list of pixels presumably 
enriched by the shower signal. For the background 
study, a similar second assigning map is created 
which contains the pixels around the signal cluster 
within a fixed margin, corresponding, in an EUSO 
like detector, to the size of two adjacent pixels. For 
each pixel of the track, the background can be esti- 
mated in two complementary ways: either by look- 
ing at the light collected by the adjacent pixels at 
the same time, or by looking at the light collected 
by the very same track pixels before and after the 
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shower signal. One obtains the observed time profile 
of the signal by adding at each time stamp the sig- 
nal from the shower pixels after subtraction of the 
background. 

4.3.2. Optics efficiency 

The focal surface may have a significant amount 
of holes which distort the detected profile as it is 
e.g. the case for the polar design (see Fig. 17b) of 
an EUSO like detector. Fig. 22 demonstrates an ex- 
ample of how the shower signal looks like on this 
type of focal surface; the corresponding distribution 
of numbers of photoelectrons as a function of time 
is shown in Fig. 24. One can notice at least three 
relatively large and rapid losses of detected photo- 
electrons (around times 13, 24 and 40 GTU) which 
correspond to holes between macrocells. The dead 
spaces between PMTs in one macrocell are smaller 
in size but also produce distortions in the profile. 

I i 
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Figure 22. An example of distribution of photoelectrons on 
the focal surface for a shower energy of 5. 10 20 eV and 
zenith angle 0=60°. Upper plot: large scale; bottom plot: 
zoomed picture, in a region where the dead space among 
PMTs modules is not negligible. 



These distortions are taken into account in the 
shower reconstruction process by building for each 
arrival time t{ an effective optics efficiency £ e ff(^)- 
The previously reconstructed 6f ov (U) and (j)tov(U) of 
the track are used to compute: 

£eff(*i) = e(x,y\9 iov (t i ),(j)f ov (t i ))A(x,y)dxdy , 



(16) 

where A(x,y) is the function returning 1 if the 
considered point with coordinates x,y is inside a 
pixel, and otherwise. e(x, y\0f ov , 0f ov ) is the optics 
point spread efficacy function on the focal surface 
at position x, y for photons incident direction deter- 
mined by 0f ov , 0f ov - The PSF and the optics trans- 
mission are taken into account to build the optics 
efficacy function. The PMT efficiencies are carried 
separately. 



4.3.3. H n 



reconstruction 



As already mentioned in section 4.2, the determi- 
nation of the detector- shower distance can be esti- 
mated in two ways: either by using the Cherenkov 
echo and the information encoded in the fluorescence 
signal profile, or, when the former is not available, 
by using fluorescence signal profile only. The latter 
case is, unfortunately, the most common (see sec- 
tion 5). In both methods, all contributions to the 
signal other than direct fluorescence and reflected 
Cherenkov lights are neglected. This allows to ob- 
tain a first estimation of the altitude of the fluores- 
cence signal maximum H max . In what follows the 
measured signal profile is corrected at each arrival 
time ti by £ e s(U) for optical loss and focal surface 
inefficiencies. 

4.3.3.1. i^max reconstruction with the Cherenkov 
peak This method needs to infer from the time pro- 
file of the signal the difference (AT) between the 
arrival time of the Cherenkov peak and that of the 
maximum of the fluorescence light t max . AT is used 
to estimate H max . 

From Fig. 23, one can find that: 



cAT 



cAT + 2R 



cos 6 
where: 

and 



2 cAT + i?(l-n max .ft) 



(17) 



R 



H — H u 



Hr 



cos(# fov (£ max )) 
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with Hq the known altitude of the reflected 
Cherenkov surface, H the telescope altitude, and 
the previously reconstructed parameters: ft the 
shower direction (with zenith angle 0) and n max the 
unit vector pointing to the telescope and defined by 

#fov(£max) and 0fov(£max). 
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Figure 23. Geometry of the Extensive Air Shower seen by 
the detector. The Cherenkov point, when available, gives the 
position on the impact point of the shower on the reflecting 
surface (ground or cloud). 



4.3.3.2. i^ m ax reconstruction without the Cherenkov 
peak The altitude of the EAS maximum can be re- 
constructed relying on the fluorescence light infor- 
mation only. The idea of the method [66] [67] is based 
on the simple observation that the time width of 
the EAS development strongly depends on air den- 
sity (and thus on the altitude) in the atmosphere at 
which EAS develops. 

A longitudinal profile of the number of shower 
electrons at given atmosphere depth X can be conve- 
niently described by an analytic function /(£) which 
depends on a dimensionless parameter £ = X ~ X1 as 



X 



n c (€) 



= e /(0 



where X\ is the depth of the atmosphere at which 
EAS starts to develop and Xq is the air radiation 
length. An example of /(£) is given in section 3.1 
for the GIL parametrization, while the considera- 
tion below is more general. The EAS profile can be 



expressed as a function of detection time t. Around 
the EAS maximum one can write 



/(c)= n|nax) (? _^ ax) 2 



f ((max)( (^max) /, , \2 
2 ^max ) 

(t tmax) 



2a 2 ' 

In this way we have obtained a gaussian distribu- 
tion of times around £ max with standard deviation a. 
In order to calculate a, we first obtain d£ / dt taking 
into account the following: 

dX _ p(h)dl _ p(h) cdt 

Xq Xq Xq 1 - 



Thus, <t is expressed by: 
X (l- 



Sl) 



a ■ 



, (18) 

P ( J fffl-)c V /-/"(Cmax) 

This general equation can be simplified further 
for the case of GIL parametrization for which 
/"((max) = — l/2£ max and the time width becomes 

-X (1 - n max • n) 



(19) 



with ( m ax only logarithmically depending on 
UHECP energy which can be easily taken into ac- 
count in an iterative procedure or fixed to a mean 
expected value. 

Because the fluorescence yield is almost constant 
between and 20 km (it varies between 4.1 and 4.7 
photons/m, see 3.2.1) and the atmospheric attenu- 
ation is a slowly varying function along the shower 
track, the shape of the shower electrons time profile 
is at first order similar to that of the total collected 
signal. This way one could measure p(H^^) and 
thus i^max of the shower maximum using observable 
informations a, n max and reconstructed ft. 

In principle, when the altitude of the Cherenkov 
scattering surface is not known, this second method 
may allow to recover it just as: 



Hn, = H- 



fluo 



Am 



4.3.4. Profile reconstruction 

The aim of this part of the reconstruction is to es- 
timate the atmospheric depth profile of the shower 
electrons from the light flux observed by the tele- 
scope. 

A functional form Fi ong for this profile is assumed, 
a possible choice is the GIL function (1), which has 
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three free parameters: E the energy of the shower, 
I max the atmospheric depth where the number of 
shower electrons reach its maximum and X\ the 
depth of the first interaction. Fi ong gives the Nf elec- 
trons of the shower as a function of the slant depth 
along the shower track. In a first estimation of E, X\ 
and X max can be set to their mean expected value. 

For each arrival time ti the expected signal at tele- 
scope Si is computed as a function of the number of 
shower electrons (unknown) to minimize the follow- 
ing x 2 function via numerical minimization meth- 
ods. 

(s>i — Si(E, X max , Xi)^j 



where Si is the measured signal (background sub- 
tracted) at time U and Gi is the uncertainty on it. 

Direct fluorescence signal observed at time ti is 
produced at a slant depth Xi in the atmosphere. 
Once a reference in altitude and Q are recon- 
structed, the EAS track in the atmosphere can be 
found unambiguously and Xi computed as an in- 
tegral along the EAS track from P^, the shower 
position at time ^, to top of atmosphere: 



Xi 



J Pi 



P(h)dl : 



this requires to know the precise density profile of 
the atmosphere for each event. 

For each time stamp the number of fluorescence 
photons Nl produced in a slant depth interval AXi 
by the Nf electrons of the shower is computed using 
the following formula: 

Nf(E,X max ,X 1 ) = 



N?(E, X max , X^a^X^, X^AXi. 



with 



AX; 



PiCt GTU 



i - ni • n 

where pi is the atmospheric density at position Pi, Yi 
is the local fluorescence yield (number of photons per 
unit of deposited energy) for atmospheric conditions 
at point Pi, di is the average energy deposit per unit 
depth per electron at shower age ^(X max ,Xi) (a 
parametrization of oti{s) can be found in [36] or set 
to its mean value of 2.2 MeV/(g/cm 2 ) ), and £gtu 
is the length of the time unit. 

The fraction Ti of photons produced at Pi that 
reach the telescope due to atmospheric attenuation 
and the distance Ri between the shower point and 



the telescope are also computed. Thus the expected 
direct fluorescence signal is expressed by: 



X max , Xi)Tie e fi(ti) 



A 



47rRj 



•2 ' 



where £ e ff(^) is the previously computed detector 
efficiency at time ti defined in (16) (see section 4.3.2) 
and A the telescope aperture. 

Although this has not yet been implemented in- 
side ESAF a more precise reconstruction process 
would require to compute the expected reflected 
Cherenkov signal, the Cherenkov scattered signal 
and the fluorescence scattered signal, but at first or- 
der the expected total signal can only be expressed 
byS/(£,X max ,Xi). 

Since the reconstructed shower track depends on 
the reference in altitude of the shower maximum, the 
fit is done iteratively until a convergence is obtained. 
Typically 2-3 iterations are enough to get a stable 
fit result. 

An example of the fit is shown in Fig. 24. 




GTU 

Figure 24. An example: Fit of temporal distribution. Points 
with error bars are observed number of p.e., blue curve - 
estimated background, red curve is the sum of fluorescence 
and background. The Cherenkov echo seen in the last bins 
of the temporal distribution was identified by the code and 
excluded from the fit. Energy and direction are the same of 
Fig. 22 



5. Simulation results 

All the simulations discussed in the present sec- 
tion are based on a standard configuration which de- 
rives from the EUSO Phase A design, i.e. a 30° FoV 
space-born telescope installed on board of the In- 
ternational Space Station at a typical orbital height 
ranging from 360 to 450 km. 
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However, depending on particular cases, the di- 
ameter of the entrance pupil is varied. In order to 
study the light flux at the pupil, independently of 
the detector, the size is fixed at 5 m. 

When the performances of an EUSO-like detec- 
tor are involved, the pupil diameter is fixed at the 
standard EUSO diameter of 2.5 m. In the trigger ef- 
ficiency studies different pupil diameters have been 
tested. 

5.1. Light flux at the pupil in clear sky conditions 

Results presented in this section concern light flux 
arriving on a 5 m pupil diameter. 

Simulation results were obtained from showers pa- 
rameterized with the GIL formula and a fluorescence 
yield based on the Kakimoto model completed with 
Bunner rays. The reference detector altitude is as- 
sumed to be 430 km and the atmosphere model is 
a US- Standard clear sky. The Earth ground is as- 
sumed to be a lambertian surface, with an albedo 
of 8%. Propagation of light is performed with the 
"Reduced" Monte-Carlo algorithm. 

Fig. 25 and 26 show an example of the simu- 
lated light at the pupil with the different light com- 
ponents: direct fluorescence, ground reflected lights 
(fluorescence and Cherenkov), Rayleigh scattered 
fluorescence and Cherenkov lights. The contribution 
of the first 4 orders of scattering are also shown. 
The signal as a function of time has been integrated 
within the whole FoV. As one can see from Fig. 25 
the direct fluorescence light contributes only to 37% 
of the total signal. This result shows the importance 
to include the scattered light contamination in the 
simulation program. The contribution of scattered 
light decreases by a factor of 2 at each scattering 
order from 35% at order 1 to only 4% at order 4. Al- 
though scattered photons are numerous, their posi- 
tions in the FoV constitute a large halo around the 
main track as shown in Fig. 26. Besides this, they 
reach the pupil later than direct signal and thus do 
not affect the first part of the profile. As a conse- 
quence, during the reconstruction process, appropri- 
ate pixels selection around the main track in space 
x-y and time will reduce the scattered light contri- 
bution as it will be explained later in section 5.1.3. 

5.1.1. Direct fluorescence light 

Fluorescence light is related to the shower devel- 
opment in atmosphere, therefore direct fluorescence 
light is the signal from which the energy and the 
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Figure 25. Number of photons reaching the pupil as a func- 
tion of time. The signal is integrated within the whole FoV. 
Top histograms refer to the last photon interaction and bot- 
tom histograms refer to the number of scattering interactions 
before reaching the pupil. The simulation has been carried 
out with 4 orders of scattering in a US-Standard clear sky 
atmosphere. The energy of the shower is 10 20 eV and the 
zenith angle is 60°. 



I max are derived. We report here its main charac- 
teristics. 

As underlined in section 3.2, inclined showers 
produce more fluorescence photons compared to 
vertical ones. Traveling in a less dense atmosphere 
these photons are moreover less scattered or ab- 
sorbed, therefore direct fluorescence light flux at 
pupil depends strongly on the shower zenith angle. 
Fig. 27 shows this dependence for showers impact- 
ing at nadir. The amount of direct fluorescence light 
increases by more than a factor 6 between vertical 
showers and horizontal ones (zenith angle above 75 
degrees) . 

At fixed zenith angle, direct fluorescence light 
reaching the telescope depends only slightly on the 
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Figure 26. Distribution of collected photons positions in the 
FoV. The simulation has been performed with 4 orders of 
scattering in a US-Standard clear sky atmosphere. The en- 
ergy of the shower is 10 20 eV and the zenith angle is 60°. 
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Figure 27. Number of direct fluorescence photons at pupil as 
a function of shower zenith angle for showers at the telescope 
nadir. Numbers of photons are normalized to the number of 
photons for a vertical shower, which is around 8000 for a 
primary energy of 10 20 eV. 

shower position in the FoV. This weak dependence 
of signal intensity within the FoV is a specific fea- 
ture of a space-based detector, due to the large 
distance between showers and telescope. Specific 
studies carried out with ESAF prove that the signal 
from showers at the edge of the FoV is reduced by 
40% compared to the nadir case. 

As already mentioned, at higher altitude, where 
the atmosphere is less dense, the shower is longer 
than at lower altitude. On the contrary, because 
the fluorescence yield is almost constant between 
and 20 km, the number of photons at the maximum 
depends weakly on the altitude. Therefore at fixed 
zenith angle the ratio between the total and the max- 
imum number of fluorescence photons is an indicator 
of the shower altitude as illustrated in Fig. 28 and 



could be used when the reflected Cherenkov signal 
cannot be observed. 
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Figure 28. Ratio of total over maximum number of fluores- 
cence photons as a function of the shower maximum altitude 
at different zenith angles. As A^ max is expressed in GTU -1 , 
the ratio of iVtot over ^max is expressed in GTU. The spread 
of the curves is due to the random first interaction point of 
UHECP used in the simulation. No correction is applied for 
the shower distance. 



5.1.2. Ground reflected Cherenkov 

Cherenkov light that is beamed around shower 
axis, travels with shower front down to Earth and 
then back to the detector for the part that is re- 
flected. It provides an extrapolation of shower axis 
to Earth, therefore the time difference between 
Cherenkov peak and fluorescence light maximum is 
an indicator of the shower altitude. 

Fig. 29 shows results obtained at different zenith 
angles. Nevertheless, this signal can only be used if 
it is sufficiently intense. 

In case of inclined showers, the Cherenkov beam 
has a longer path to go down to the ground and 
is therefore more attenuated than in the case of 
vertical showers. Moreover Cherenkov photons are 
more dispersed around the shower impact due to 
the photon angular distribution. As a consequence 
the Cherenkov peak amplitude decreases when the 
zenith angle increases (Fig. 30). Being quite con- 
stant up to 40° the signal is reduced by a factor of 
30 at 75°. 

In most cases Earth surface cannot be modeled 
by a Lambertian surface and specular reflection has 
to be used. 

Taking into account showers with an impact at 
nadir, in the case of a specular reflection on oceans 
with a wind speed of 2 m/s (respectively 9 m/s), 
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Figure 29. Time difference between ground reflected 
Cherenkov and direct fluorescence peaks as a function of the 
altitude of the shower maximum at different zenith angles. 
The spread of the curves is due to the random first interac- 
tion point of UHECP used in the simulation. No correction 
is applied for the shower distance. 
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Figure 30. Amplitude of the ground reflected Cherenkov peak 
as a function of shower zenith angle. Numbers of photons are 
normalized to the number of photons for a vertical shower, 
which is around 2900 for a primary energy of 10 20 eV. 
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Figure 31. Amplitude of the ground reflected Cherenkov 
peak as a function of the shower impact position in the FoV 
for showers with a zenith angle of 20° in case of specular 
reflection over ocean with a wind speed of 2 m-s — 1 . The total 
albedo is 8%. 

have a reflected Cherenkov light lower than the one 
expected in the Lambert ian case. Furthermore, for 
showers with zenith angle greater than 30° , the spec- 
ular peak does not point toward the telescope, re- 
gardless of the shower impact position in the FoV. 
Over the oceans the hypothesis of a Earth Lam- 
bertian surface seems to be too optimistic and the 
anisotropy of the specular reflection will be responsi- 
ble of a dependence of the reflected Cherenkov light 
with the position in the FoV. 

We conclude, therefore, that in case of specular 
reflection on Earth, the Cherenkov light of most 
of the showers will not be detected. Simulation re- 
sults have been produced with a constant albedo 
of a few percents. In case of specular reflection, 
the Cherenkov peak could be enhanced only if, by 
chance, the shower angle and its impact position in 
the FoV allows the light to be reflected towards the 
telescope. 



one can deduce from Fig. 9 that events with a 
zenith angle greater than 20° (respectively 30°) will 
give a weaker Cherenkov peak than in Lambertian 
case. On the contrary, for smaller zenith angles, the 
Cherenkov peak should be greatly enhanced. 

The anisotropy of the specular reflection has 
consequences for the dependence of the reflected 
Cherenkov light on the shower position in the FoV. 

As an example, Fig. 31 shows the amplitude of 
the ground reflected Cherenkov peak as a function 
of the position of the shower impact in the FoV for 
showers with 20° zenith angle. For a wind speed of 
2 m/s, 76% of the showers that impact in the FoV 



5.1.3. Scattered light 

In this section, the term "scattered light" repre- 
sents all photons that have been scattered by Earth 
surface (assumed lambertian) or by atmosphere in 
clear sky conditions. The signal component previ- 
ously named "ground reflected Cherenkov" is not 
included in this scattered light component, but is 
considered as a component of the basic signal. As 
already said, scattered photons are slightly time 
delayed compared to the direct signal and spread 
around the main track in the FoV. Therefore they 
can partly mix with the basic signal as a noise. 

Fig. 32 shows time profiles of photons reaching the 
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Figure 32. Number of photons reaching the pupil as a func- 
tion of time at 4 zenith angles. The signal is integrated in 
the whole FoV. In gray are the basic signals: direct fluores- 
cence and ground reflected Cherenkov. The single scattered 
signal is superimposed in light red and the multi-scattered 
signal in dark blue. 



pupil at different zenith angles. The plotted signal is 
integrated in the whole FoV. Although the scattered 
light increases with shower zenith angle (by a factor 
of 5 for single scattered photons between 0° and 80° 
and by a factor of 4 for multi-scattered photons), its 
contribution is of the same order of magnitude as the 
direct fluorescence light. However part of the scat- 
tered signal is not correlated with the track in the 
FoV because most of the photons come from a last 
scattering position (in the atmosphere) far from the 
shower axis. Thus most of the scattered photons will 
be rejected in the same way as background photons. 
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Figure 33. Integrated single scattered photon number as a 
function of the distance to the shower axis for a vertical 
shower. Fluorescence and Cherenkov photons are compared. 
Distributions are normalized to 1. 
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Figure 34. Contributions of selected scattered photons from 
different origins related to direct fluorescence light as a func- 
tion of shower zenith angle in clear sky conditions. Selected 
scattered photons are those correlated in space and time 
with the direct fluorescence signal. 

As an example Fig. 33 shows the lateral distribu- 
tion of the last scattering position in the atmosphere 
for vertical showers. As one can see, due to the fact 
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Figure 35. Contributions of scattered photons to the integral 
of the signal iVtot, the integral of the first half part of the 
signal A/]/2 and the maximum of the signal iV max . 

that Cherenkov light is beamed at emission, 90% of 
single scattered Cherenkov photons have their last 
scattering position in the atmosphere located less 
than 0.4 km from the shower axis, in contrast to only 
10% of single scattered fluorescence photons. Sim- 
ulations performed to study multiple scattering ef- 
fects prove that only 5% of the photons scatter more 
than once before reaching the detector. 

If not taken into account during the reconstruc- 
tion process, the scattered light component would 
entail systematic errors on the shower parameters. 
The amount of light is increased and the detected 
shower profile is distorted. 

To estimate precisely the contribution of scattered 
light to the detected signal, we have selected the 
scattered photons, at the photo-electron level, using 
a simple detector simulation. It uses EUSO config- 
uration (pupil diameter of 2.5 m), applying a global 
efficiency coefficient (5%) to convert photons into 
photo-electrons. The FoV being pixelized and time 
profile being discretized in GTU, one selects at a 
given GTU the scattered signal standing in pixels 
containing the direct fluorescence signal. 

It appears that the contribution of scattered pho- 
tons represents less than 30% of the direct fluores- 
cence light. This is shown in Fig. 34 where contri- 
butions of scattered photons from different origins 
are represented as a function of the shower zenith 
angle. The main contribution is that of single scat- 
tered Cherenkov photons that represents slightly 
more than 20% of the direct fluorescence photons, in 
contrast to only less than 4% for single scattered flu- 
orescence photons and 1% for multi-scattered pho- 
tons. 

It has to be noticed that presence of scattered pho- 



tons has consequences not only on the total number 
of photons reaching the pupil but also on the shape 
of the reconstructed profile. This can be illustrated 
by comparison of scattered photon contributions to 
the integral of the signal, to the integral of the first 
half part of the signal and to the maximum of the 
signal as shown in Fig. 35. 



5.2. Triggering efficiency and energy threshold 

All the simulations discussed in the present sec- 
tion are based on a standard configuration which 
derives from the EUSO Phase A design at a typical 
orbital height H = 400 km. 

From this basic configuration, four different de- 
tectors were implemented, identical in all parts but 
for the diameter of the entrance pupil: 2.5 m (EUSO 
Phase A configuration), 4 m, 8 m and 12 m. In this 
way ESAF is able to evaluate the performances of 
EUSO and of possible larger detectors. 

A sample of 3-10 5 simulated EAS with energy 
in the range 10 19 eV < E < 5-10 20 eV (flat energy 
distribution in a logarithmic scale and angular dis- 
tribution uniform in solid angle) was used for all the 
configurations. 

The showers were generated by GIL (see sec- 
tion 3.1) in clear sky and US- Standard Atmosphere 
(no clouds). A parameterized optics with a 2D gaus- 
sian PSF of 5 mm RMS, a GTU of t GTV = 2.5 fis 
and an overall photo-detector efficiency of 10% were 
used according to the EUSO Phase A design. 

Because of the limits imposed by the available 
telemetry and computing resources in the EUSO 
project, only the contiguity tracking triggers that 
satisfy ^ftr ^ 100 Hz were considered useful. Fu- 
ture detectors may have less stringent requirements. 

As shown in Fig. 36, if the condition z^ftr ^ 100 Hz 
is not satisfied, the efficiency curve flattens at low 
energies to a value that depends on [L; k) (see sec- 
tion 3.4.4 for the definition of [L; fc]), else flattens to 
zero at low energies. Furthermore, when the value 
of L is raised to a fixed /c, the efficiency curve is 
shifted toward higher energies (without changes of 
shape), i.e. the full efficiency is obtained at a larger 
value of E. 

The dependency of the trigger efficiency on the 
zenith angle is shown in Fig. 37. The efficiency 
curve is shifted toward higher energies for < 30° 
and toward lower energies for > 60°. 

To compare the detectors on the same grounds, 
the continuity trigger pattern must satisfy ^ftr ^ 
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Figure 36. Trigger efficiency curves for various trigger config- 
urations (2.5 m detector): (a) FTR above 100 Hz; (b) FTR 
below 100 Hz (c) Efficiency with fixed L = 6 and 3 < k < 6. 

100 Hz for all detectors. [L;k] = [3; 4] was the se- 
lected trigger for the further studies being the less 
restrictive pattern satisfying this additional condi- 
tion. 

The results are shown in Fig. 38. As expected the 
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Figure 37. Efficiency of the [3;4] and [4;4] trigger patterns on 
three ranges in 0. The relative weight of the three samples 
are (a) ~ 25%, (b) ~ 50%, (c) ~ 25%. The red line at 
10 20 eV is drawn as a reference. 

curves are shifted toward lower energies when in- 
creasing the entrance pupil diameter. The energy 
threshold of the experiment scales as the inverse of 
the diameter or, equivalently, as the inverse of the 
square root of the photo-detection efficiency. 

Simulations with ESAF show that an EUSO 
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Figure 38. Trigger efficiency curves for various sizes of the 
detector using the pattern [3;4]. 

like detector has full detection efficiency around 
2. 10 20 eV and 50% efficiency at 8. 10 19 eV. Better 
efficiency can be obtained only by enlarging the 
detector size or by using better photo-detection 
devices than commercially available MAPMTs. We 
recall, however, that these results are obtained in 
ideal conditions, not taking into account variable 
weather and dead pixels in the focal surface. 

This result is really very weakly dependent on the 
specific triggering algorithm. The lack of efficiency 
is related to lack of signal at low energy. 

It is worth recalling that these conclusions apply 
to photo-detector similar to MAPMTs, that have 
a global photo-detection efficiency (quantum effi- 
ciency x PMT practical factor x geometric factors) 
of about 10%. To our knowledge, better devices do 
not exist yet, but new photocathodes with much im- 
proved quantum efficiency are already available and 
might yield more efficient MAPMTs in the near fu- 
ture. This of course may have a significant impact 
on efficiency. However, higher photo-detection effi- 
ciency also means higher nightglow background, so 
the gain would not be linear. To be more precise, 
the gain in photon collection efficiency scales as the 
pupil area in Fig. 38. For instance, referring to a 2.5 
m diameter pupil (right-most curve in Fig. 38), an 
increase of efficiency of a factor 2.56 (significant im- 
provement!) corresponds to 4 m diameter (second 
curve from right). 

5.3. Angular resolution 

Fig. 39 shows the angular resolution concerning 
the identification of the shower axis direction in 
space in presence of typical random background 



of 500 ph m _2 ns _1 sr _1 (see section 3.5) 10 ; in the 
figure, ^68 is the angle in space between the sim- 
ulated and reconstructed directions which defines 
the cone containing 68% of the events, i.e. one stan- 
dard deviation in angle; 0\ oc is the zenith angle of 
the incoming shower calculated at the ground im- 
pact point. Standard EUSO parameters are used 
for this simulation. The showers are simulated with 
UNISIM in clear sky conditions neglecting the mul- 
tiple scattering effects. 
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Figure 39. Angular resolution. ^68 is the angle in space 
between the simulated and reconstructed directions which 
defines the cone containing 68 % of the events. 

The curves labeled NE1 refers to exact formu- 
las and numerical methods, while AA1 refers to ap- 
proximated formulas and analytical methods (see 
section 4.2). As expected the direction reconstruc- 
tion is far better for very inclined showers, since 
for such inclinations the track length is greater and 
consequently the total number of detected photons 
and the number of hit pixels are higher (see sec- 
tion 5.1.1). 

At 10 20 eV the angular resolution is close to 10° 
for zenith angles smaller than 40°, it goes down to a 
few degrees in the angular range 50°-70° and finally 
approaches about 1° for very inclined showers with 
zenith angle greater than 80°. The angular resolu- 
tion reaches 1° already at 70° zenith angle for 2 10 20 
eV and at 60° zenith angle for 5 10 20 eV. 

In a conservative way, the detection of Cherenkov 
light is not taken into account here and the detector 
shower distance is set in the reconstruction at the 
expected average value. Shower to shower fluctua- 



10 We recall here that the typical night glow background in 
EUSO is 500 ph m^ns -1 sr _1 , reaching 1000 with 1/4 of 
moon. With more than 1/4 moon, the observation of EAS 
is not possible. 
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tions yields approximately 0.5 km rms error on the 
altitude. 

The relation between the uncertainty in the es- 
timation of the shower maximum height (and so of 
the detector-shower distance) and the uncertainty 
in the reconstructed direction has been studied by 
Monte- Carlo simulation and analytical calculation. 
In average we find an uncertainty of ~0.2° in the 
reconstructed zenith angle for every km of error in 
the altitude of the EAS maximum. 



5.4. Energy and X max reconstruction performances 

As previously seen, the event reconstruction starts 
with the determination of the shower track. Then 
a functional form for the shower longitudinal devel- 
opment is assumed and the expected collected light 
is calculated. The energy and X max of the EAS are 
deduced by varying the parameters of the longitu- 
dinal shower function until an agreement with the 
measured signal is obtained. 

Several sources of uncertainties can be identified 
at each step of this reconstruction process. ESAF 
is a useful tool to evaluate most of them. Different 
samples of proton showers have been generated. In 
each sample, either one shower parameter (the en- 
ergy, the zenith angle, the shower maximum alti- 
tude, the ground impact position in the FoV) or the 
atmosphere profile, or one telescope characteristics 
is varied. Possible uncertainty sources can be stud- 
ied using each of these simulated samples. 

5.4.1. Uncertainties on the shower track altitude 

Beside the angular resolution described in sec- 
tion 5.3, the precision on the shower track geome- 
try is related to the knowledge of the shower track 
altitude. As explained, two different methods may 
be applied to determine i^max, either using the re- 
flected Cherenkov signal, or when the information 
can not be retrieved, using the shower track length 
only. In both cases, the shower track zenith angle 
must be known to infer H max . 

At the first order, according to geometrical con- 
siderations only, the impact on i^ max of an uncer- 
tainty on is estimated. Whereas an uncertainty 
SO = +1° induces an underestimation of i^ max when 
the Cherenkov information is used (relatively in- 
creasing with but less than 2% for shower with 
zenith angle < 60°), it induces an overestimation 
of i^max when the Cherenkov information is not used 



(relatively decreasing with but less than 2% for 
shower with zenith angle > 20°). 

The results are obtained simulating about 12600 
events (of which 11200 reconstructed) in clear 
sky (US-Standard atmosphere) with L0WTRAN and 
UNISIM and a telescope aperture of 2.5 m (para- 
metric optics). 

ESAF results show that in clear sky and US- 
Standard Atmosphere, i^max can be reconstructed 
with an error of 100-300 m for events with both 
fluorescence and Cherenkov echo components for 
10° < < 60°, and with an error of 400-800 m when 
the Cherenkov echo is either missing or just ignored 
for 10° < < 90°. The absolute error grows with 
for both methods, while the relative error is rising 
for the method based on the Cherenkov echo (from 
5% to 8%) and decreases for the method based on 
the time width of the p.e. distribution (from 22% 
to 5%). Therefore these methods can be considered 
as complementary. 



5.4.2. Effects of geometric uncertainties 

Uncertainties on the shower geometry affect the 
factors used in the expected light profile calcula- 
tion. At each shower point P^, the reconstruction 
process involves the fluorescence yield the light 
transmission T^, and a solid angle factor (AirRf) -1 . 
ESAF allows to evaluate uncertainties induced on 
them related to shower geometry. A ~1 km shift on 
altitude affects the solid angle computation by less 
than 1%. Despite the fact that the light transmis- 
sion is the more affected factor by an uncertain al- 
titude, an uncertainty on the reconstructed altitude 
of ~1 km leads to less than 7% uncertainty on the di- 
rect transmission in clear sky. The fluorescence yield 
dependence with altitude is almost canceled by the 
dependence of the ionization energy deposit per me- 
ter of track length, thus an uncertainty of ~1 km on 
altitude affects the fluorescence yield per meter by 
less than 1%. 

Furthermore, the uncertainty on the shower ge- 
ometry affects the slant depth calculation. The inte- 
gral of the reconstructed profile is thus affected by 
an uncertainty of 10% for every km of error in alti- 
tude, and X max by around 100 g/cm 2 . In the same 
way an uncertainty on shower zenith angle leads to 
an uncertainty on X max of around —20 g/cm 2 for 
SO = +1° but to only 0.2% for SO = +1° on the 
integral of the shower profile. 
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5.4.3. Effects of atmospheric uncertainties 

ESAF allows also to evaluate uncertainties on Y i: 
Ti and slant depth calculation related to atmosphere 
knowledge. During simulation runs, atmosphere pa- 
rameters are varied according to all available MSISE 
profiles, and in the reconstruction process, only US 
Standard atmosphere model is considered. These 
studies show that the uncertainty on transmission 
and on the fluorescence yield per meter stays below 
2% (respectively 2.5%). They also demonstrate that 
uncertainties induced by slant depth calculation lead 
to a shift on X max from —30 g/cm 2 to 25 g/cm 2 , 
and to an uncertainty on the reconstructed profile 
integral around 3% which could reach 4% for shower 
detected close to the equator. 



also due to the increasing number of collected pho- 
tons as 6 increases. 

As already said, an uncertainty on the shower al- 
titude induces an important uncertainty on the re- 
constructed X max , that is why results on X max reso- 
lution are presented depending on the method used 
for the first altitude estimation. For events with both 
fluorescence and Cherenkov signals X max can be re- 
constructed with resolution better than 40 g/cm 2 
for < 60°. At larger angles the Cherenkov signal 
becomes suppressed, and the method based on the 
fluorescence profile must be used. In this case the 
resolution on X max is degraded, but it is still in the 
range between 50 g/cm 2 and 100 g/cm 2 , depending 
on the shower inclination. 



5.5. Energy and X max resolution 

Apart from the uncertainties on the atmospheric 
profile, different contributions described above are 
taken into account to perform preliminary studies of 
the relative resolution of energy a/E for an EUSO 
detector with 2.5 m entrance pupil diameter. Show- 
ers are simulated with UNISIM, standard EUSO pa- 
rameters are used for the detector and simulations 
have been performed in clear sky conditions with 
a US-Standard atmosphere, neglecting the multiple 
scattering effects. Results are displayed in Fig. 40 as 
a function of UHECP energy in various bins of zenith 
angle for proton initiated EAS. An energy resolu- 
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Figure 40. Relative energy resolution a/E as a function of 
energy in various bins of zenith angle for proton initiated 
EAS. 

tion between 20-25% is expected at E = 10 20 eV. It 
worsens at lower energy due to a low signal statistics 
and improves at higher energies. There is a strong 
dependence on the zenith angle of developing EAS, 
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Figure 41. The plot shows the error (in g/cm 2 ) with which 
I max is reconstructed when the Cherenkov pulse is available 
(lower curve, blue points) and when only fluorescence profile 
is used (upper curve, red points). 

5.5.1. Effect of scattering on E and X max resolution 
Simulation results show that photon scattering 
distorts the fluorescence signal. Studies described 
in section 5.1.3 are used to evaluate systematics re- 
lated to the fluorescence signal contamination. In 
case of clear sky conditions, the total number of 
photons at the entrance pupil is increased by 10% 
to 25%, depending on the shower zenith angle (see 
Fig. 34). The contamination is mainly due to the 
scattered Cherenkov light. If this contribution is not 
subtracted to the signal, it leads to an uncertainty 
on E in the same range, and to a shift of the max- 
imum of the profile (in time); X max is then overes- 
timated by around 100 g/cm 2 . Scattering of fluo- 
rescence photons as well as multi-scattering (more 
than once) of the Cherenkov photons have weaker 
impacts on the energy (less than 5%) and X max 
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(10 g/cm 2 ) resolutions. The subtraction of the dif- 
fused Cherenkov light is possible, as it was already 
done in HiRes experiment (see [68]). 



5.6. Effect of clouds on the detected signal 

5.6.1. Reflected peak 

Because clouds and aerosols are responsible for 
Mie scattering, an adequate simulation including 
multi-scattering must be carried out. As examples 
Fig. 42 and 44 show simulation results at the pupil 
entrance of the same shower as in section 5.1 but 
with presence of rural aerosols or clouds. As one 
can see, contributions from scattering orders greater 
than 4 cannot be neglected as in the clear sky case. 
Therefore we carried out the simulations taking into 
account up to 8 scattering orders. 

In presence of aerosols, fluorescence photons be- 
ing emitted mostly above aerosols layers, direct flu- 
orescence light is only weakly affected. On the con- 
trary, ground reflected light is strongly attenuated 
since it has to travel twice through the dense layers 
of atmosphere before being detected. Aerosols be- 
ing mostly a scattering medium, a multi-scattering 
peak dominated by the contribution of Mie scat- 
tered Cherenkov photons appears in the total pho- 
ton time distribution (Fig. 42). Its maximum is be- 
tween time corresponding to ground impact and that 
corresponding to aerosol top layer impact. However 
most of these photons are not contiguous with the 
shower track. To estimate the shape of the multi- 
scattering peak likely to be detected, an appropriate 
selection of the scattered signal component at the 
photo-electron level is applied, using a simple detec- 
tor simulation (already mentioned in section 5.1.3). 
Only the scattered signal which is contiguous in 
space and time with the shower track is selected, 
keeping in mind that selecting too many pixels re- 
duces the signal to noise ratio. 

This selection has been applied to the considered 
event and the resulting signal is shown on Fig. 43. 
As one can see the amplitude of the multi-scattering 
peak is weak compared to statistical fluctuations. 
Let us notice moreover that the peak maximum af- 
ter pixel selection in space and time still does not 
correspond to the impact on the aerosol top layer. 
In a general way, its position in time should depend 
on the aerosols properties. Therefore the use of the 
multi-scattering peak (if detected) instead of ground 
reflected Cherenkov appears in the reconstruction 
procedure to be more difficult than expected. 
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Figure 42. Number of photons reaching the pupil as a func- 
tion of time. Signal integrated in the whole FoV. Top his- 
tograms refer to the last photon interaction and bottom his- 
tograms refer to the number of scattering interactions before 
reaching the pupil. The simulation has been carried out with 
8 orders of scattering in a US-Standard atmosphere includ- 
ing a rural aerosol model of visibility 5 km. The energy of 
the shower is 10 20 eV and the zenith angle is 60°. 



Contrary to aerosols, clouds can be found at 
higher altitudes, and thus distort strongly direct 
fluorescence light depending on their altitude, op- 
tical depth and the shower zenith angle. However, 
in some cases of low altitude clouds, when major 
part of fluorescence light (including the maximum 
of development) is emitted above clouds, direct flu- 
orescence light is only weakly distorted. Neverthe- 
less reflected Cherenkov is strongly attenuated and 
scattered photons are added to signal. As in case 
of aerosols, a multi-scattering peak dominated by 
Mie scattered Cherenkov photons appears in total 
photon time distribution (Fig. 44). Mie scattered 
photons are also spread in FoV therefore the ampli- 
tude of the multi-scattering peak is greatly reduced 
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Figure 43. Mean photoelectrons number as a function of time 
after selection of the shower track contiguous pixels. The 
energy of the shower is 10 20 eV and the zenith angle is 60°. 
The atmosphere model is US-Standard with rural aerosols 
of visibility 5 km. The vertical bars indicate respectively the 
impact on aerosols top layer and the impact on ground. 

and low compared to statistical fluctuations when 
pixel selection is applied (Fig. 45). We have studied 
several configurations as a function of shower zenith 
angle, cloud top altitude, cloud thickness and cloud 
optical depth. In some particular situations, for 
dense clouds located at low altitude and showers 
weakly inclined (< 40°), the multi-scattering peak 
may be clearly visible. Two peaks may even appear 
in very specific situations, the first related to the 
impact on the cloud layer, the second related to 
the impact on ground. In a general way, the peak 
shape and its position in time depend strongly on 
clouds properties. Therefore, in the reconstruction 
procedure, the use of the multi-scattering peak (if 
detected) instead of ground reflected Cherenkov 
appears to be more difficult than expected. 

An on board atmospheric sounding device associ- 
ated with the telescope is thus required in order to 
interpret the signal distortions due to aerosols and 
cloud layers. 

5.6.2. Effect of clouds on acceptance 

Presence of clouds and aerosols in the FoV makes 
signal more complex to exploit and modifies also 
detector acceptance. Despite the fact that multi- 
scattering until order 8 is needed in cloud simula- 
tions and because the "Reduced" Monte-Carlo al- 
gorithm for radiative transfer is time-consuming, a 
compromise using the "Bunch" algorithm has been 
done for preliminary studies concerning acceptance 
with clouds. The main effect of clouds on the accep- 
tance comes from the attenuation of the detected 




Time (GTU) 

Figure 44. Number of photons reaching the pupil as a func- 
tion of time. Top histograms refer to the last photon interac- 
tion and bottom histograms refer to the number of scattering 
interactions before reaching the pupil. The simulation has 
been carried out with 8 orders of scattering in a US-Stan- 
dard atmosphere including an homogeneous cumulus cloud 
model, characterized by a top altitude of 3 km, a vertical 
extension of 1 km and an optical depth of 2. The shower 
energy is 10 20 eV and the zenith angle is 60°. 

fluorescence light when a cloud layer stands between 
a shower and the detector. For this study, the scat- 
tering component has not been calculated, only the 
attenuation effect has been considered, using a sim- 
plified version of the "Bunch" algorithm. 

Two samples of 2-10 4 proton EAS with flat energy 
distribution (in a logarithmic scale between 10 19 eV 
and 10 21 eV), isotropic angular distribution at top of 
atmosphere and a maximum of shower development 
in FoV have been used. Showers were generated with 
GIL parametrization in an atmosphere distribution 
obtained from MSISE model combined with the In- 
ternational Space Station Trajectory data. TOVS 
distribution and description of clouds (probability of 
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Figure 45. Mean photoelectron number as a function of time 
after selection of shower track contiguous pixels. The shower 
energy is 10 20 eV and the zenith angle is 60 degrees. The 
atmosphere model is US-Standard with an homogeneous cu- 
mulus cloud model, characterized by a top altitude of 3 km, 
a vertical extension of 1 km and an optical depth of 2. The 
vertical bars indicate respectively the impact on cloud layer 
and the impact on ground. 
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Figure 46. Trigger efficiency curves in a modeled atmosphere 
with and without clouds. 

presence, opacity and altitude) have been included 
in one of the two samples. The Earth surface is con- 
sidered as a lambertian one, with an albedo of 5%. 
The fluorescence yield is based on the Kakimoto 
model completed with Bunner rays. The detector is 
derived from the EUSO Phase A design and the trig- 
ger efficiency is computed using the pattern [4;4]. 

In Fig. 46 the effect of clouds on the trigger effi- 
ciency is illustrated: it decreases when presence of 
clouds is taken into account. At the highest energies, 
between 10 20 eV and 10 21 eV, the mean efficiency is 
only 80% when clouds are taken into account instead 
of 96% in clear sky conditions. 
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Figure 47. Acceptance for triggering showers that have their 
maximum within the FoV as a function of the zenith angle 
in clear sky conditions (top) and cloudy conditions (bottom). 
The plot also shows the contribution of events with (in blue) 
and without (in red) Cherenkov peak. The shower energy 
spectrum is flat (in a logaritmic scale) between 10 20 eV and 
10 21 eV. To guide the eyes, the curves giving an estimate of 
the geometric acceptance (showers having their maximum in 
the FoV) are also superimposed. Percentage are related to 
the simulated geometric acceptance whose integrated value 
is 633 000 km 2 sr (the estimated one gives a 4% lower value 
than the simulated one). 



In Fig. 47 we also show the trigger acceptance as 
a function of the shower zenith angle in clear sky 
conditions and in the presence of an atmosphere 
with clouds. In clear sky conditions the data nicely 
agree with the expectation that is at first order pro- 
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portional to sin(20). The contribution to the total 
number of events of near horizontal proton initiated 
showers is negligibly small; however, the possibility 
to reconstruct and identify such showers might pos- 
sibly open a new window for the detection of high 
energy neutrinos. We have not considered in this pa- 
per the problems related to neutrino detection. 

The main effect of clouds, as already said, is that 
the global trigger acceptance is decreased from 96% 
to 80%, but also that events having a detectable 
Cherenkov peak is reduced from 79% down to 36%. 
To obtain this last result we have defined a "de- 
tectable Cherenkov peak" as a peak with maximum 
amplitude greater than 4 standard deviations above 
background. 

In case of clouds presence, a triggered event may 
result in a strongly distorted signal. To estimate the 
rate of showers with a maximum of development 
standing above cloud layer, if any, a condition on 
shower age has been applied to simulated events: 
only events for which the portion of the longitudinal 
profile between s = 0.9 and s = 1.1 develops above 
potential clouds are kept. As shown in Fig. 48, this 
condition reduces the efficiency by a factor of 2. This 
result is a mean value assuming a realistic cloud 
coverage. 
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Table 1 

Events collected in 3 years of data taking for various en- 
ergy intervals by an EUSO like detector and by a detector 
with a pupil 4 times larger with the assumption of a cosmic 
ray spectrum as measured by HiRes collaboration [69]. The 
first two columns report the triggering events. The second 
two report the events that have their maximum well above 
a cloud layer with average world cloud coverage. The condi- 
tions and the assumption made to obtain these numbers are 
explained in text. 
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Table 2 

Same as in Table 1 but numbers are obtained for a cosmic 
ray spectrum modeled by a single power law with spectral 
index equal to 2.75 

flux of 3. 10 -33 m -2 s _1 sr _1 eV -1 at the energy 
of 10 19 eV. The behaviour of this spectrum was re- 
cently confirmed by the Pierre Auger collaboration. 
The number of generated protons takes into account 
a duty cycle of 18% and the acceptance of 633000 
km 2 sr. Simulation conditions are the same as pre- 
viously mentioned. Results are presented in Table 1. 
In Table 2 we also report the same data as in Table 
1 but using a different energy spectrum, namely the 
one that was suggested by Agasa data at the time 
when the EUSO telescope was studied. 



Figure 48. Trigger efficiency in a MSISE modeled atmosphere 
with TOVS clouds. Superimposed in red: triggered events 
which have part of their longitudinal profile (between s = 0.9 
and s = 1.1) located above clouds. 

To estimate the event rates for 3 years of data tak- 
ing of an EUSO-like experiment, 13316 protons have 
been simulated between 4.10 19 and 10 21 eV. The en- 
ergy spectrum is given by [69] and normalized to a 



6. Conclusions 

In order to study the scientific potential of an ob- 
servatory in orbit around the Earth, a full simulation 
of the detection of EAS resulting of UHECP inter- 
action in atmosphere has been performed. Although 
this work has been done within the framework of 
the ESA EUSO project, the results are quite gen- 
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eral and apply to a large class of possible telescopes. 
The code also implements a complete reconstruction 
chain of direction, energy and shower X max . 

The simulation of the radiative transfer in the at- 
mosphere up to the telescope emphasizes the im- 
portance of multiple scattering in atmosphere and 
Earth. The contamination of the direct fluorescence 
signal by the scattered light is enhanced, particu- 
larly in presence of clouds or aerosols: in these cases 
the longitudinal profile is distorted. The behaviour 
of the reflected Cherenkov light, according to the 
ground and atmospheric parameters, has been stud- 
ied. The simulation results show that the effects of 
the atmosphere on the EAS light signal have to be 
considered during the reconstruction procedure. 

The studies reported in this paper show that the 
detection energy threshold (defined as the energy 
above which detection efficiency exceeds 50%) can- 
not be significantly less that 10 20 eV, unless the aper- 
ture is much larger than the one originally foreseen 
for EUSO (2.5 m in diameter). However, it is shown 
that a large EUSO-like telescope of the order of 8 to 
12 m in diameter might perform very well around 
the GZK region. In fact, while a triggering threshold 
of 2. 10 19 eV can be obtained only with a pupil di- 
ameter of at least 8 m with current photo-detection 
technology. 

The measurement of the arrival direction of the 
UHECP depends mostly on the focal surface gran- 
ularity and on the electronics time resolution, but 
with a small detector like EUSO the angular resolu- 
tion is also limited by lack of signal. Below 10 20 eV 
the angular resolution is always worse than 1° even 
for very inclined showers and is about 4° on average 
at 10 20 eV between 30° and 70°. 

Significant improvements are possible only with 
a detector larger than EUSO (larger signal with 
smaller fluctuations) and possibly with a high gran- 
ularity focal surface. 

The ESAF code, the tool that we have developed 
to make this study, is available upon request and is 
suitable for future developments in this field. 
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